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ABSTRACT 

Recent  studies  have  shown  that  various  inaccuracies  occur  when 
a  single -degree -of -freedom  gyroscope  is  subjected  to  an  angular  vibra- 
tory environment.  Under  such  conditions  the  performance  of  the  gyro- 
scope has  been  analyzed  and  the  nature  of  the  inaccuracies  investigated. 
Errors  due  to  "coning",  "dynamic  unbalance",  and  "sculling"  have  been 
included  in  the  investigation. 

Values  of  the  inaccuracies  representative  of  those  expected  from 
tests  using  the  precision  angular  vibrator  developed  by  the  Instrumentation 
Laboratory  of  the  Massachusetts  Institute  of  Technology  are  presented. 

By  reference  to  a  dynamic  model  the  equation  of  motion  was  ob- 
tained.   A  suitable  solution  to  the  equation  was  determined  in  series 
form  for  the  case  of  sinusoidal  inputs  simultaneously  applied  about  two 
axes,    No  restrictions  are  placed  on  the  frequency  of  the  applied  mo- 
tions.   The  terms  in  the  series  solution  reveal  the  nature  of  the  angular 
displacement  of  the  float  with  respect  to  inertial  space  about  the  output 
axis  for  two  cases  :    (1)  zero  elastic  restraint,  (2)  finite  elastic  restraint. 

It  is  strongly  recommended  that  correlation  of  actual  experimen- 
tal results  with  the  developed  theory  be  obt  ained  prior  to  further 
theoretical  considerations. 


Thesis  Supervisor:    Dr.  Walter  Wrigley 

Title:    Professor  of  Instrumentation  and 
Astronautics 
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OBJECT 

The  object  of  this  thesis  is  to  analyze  the  performance 
of  a  floated  single -degree -of -freedom  integrating  gyroscope 
subjected  to  angular  vibrations  in  order  to  determine  the 
inaccuracies  that  may  result  from  such  an  environment. 


CHAPTER  1 


INTRODUCTION 


The  age  of  the  golden  sixties  is  upon  us.    Ushered  in  with  this 
decade  is  an  era  of  accomplishment  which  has  raised  the  importance 
of  the  gyroscope  to  an  all-time  high.    This  era  of  great  advancement 
includes  those  fields  of  Science  and  Engineering  which  are  brought  to 
the  fore  when  one  thinks  of  space.    Because  of  space,  these  fields 
have  been  searched  diligently  for  methods  of  navigation  and  control 
and  from  this  searching  the  gyroscope  has  emerged,  holding  a  posi- 
tion of  extreme  prominence. 

With  this  relatively  new  stature,  the  gyroscope  has  naturally 
been  investigated  in  efforts  to  improve  its  performance  as  well,  as 
to  predict  what  may  be  expected  in  the  future.    The  floated  single- 
degree -of -freedom  integrating  gyroscope,  an  important  element  in 
most  space  navigation  and  guidance  systems,  is  the  subject  of  this 
analysis.    The  SDFI  gyroscope  has  been  interrogated  from  the 
standpoint  of  being  subjected  to  varying  vibratory  inputs  which  may 
introduce  error.    These  vibrations  may  be  of  very  low  frequency 
due  to  the  bending  modes  of  long  missile  bodies  or  they  may  be  of 
high  frequencies  associated  with  aerodynamic  flutters  and  buzzing. 
Through  the  theoretical  results,  it  is  hoped  that  undesirable  effects 


caused  by  such  vibratory  motions  will  be  recognized  and  later  minimized 
when  actual  verifying  data  is  available. 

The  floated  single -degree -of -freedom  integrating  gyroscope  is 
an  instrument  capable  of  sensing  angular  velocities  to  a  high  degree 
of  accuracy.    When  utilized  in  inertial  guidance  systems,  the  gyro- 
scope must  have  considerable  stability  and  extreme  sensitivity.    The 
vibratory  environment  in  which  the  gyro  functions  may  introduce  un- 
wanted angular  outputs.    This  degradation  of  desired  performance  is 
the  result  of  both  linear  and  angular  oscillations.    An  analysis  of  the 
undesired  effects  is  made  in  order  to  determine  the  nature  and  magni- 
tudes of  the  errors  involved. 

This  investigation  was  conducted  with  the  view  that  the  Massachu- 
setts Institute  of  Technology  Instrumentation  Laboratory  will  conduct 
tests  on  a  precision  angular  vibrator.    A  previous  investigation  has 
disclosed  what  may  be  expected  if  two  equal  frequencies  are  simul- 
taneously applied  to  the  output  and  spin  reference  axes  of  a  typical 
single -degree -of -freedom  gyroscope.    The  present  analysis  is  concerned 
with  any  two  frequencies  simultaneously  applied  to  the  output  and  spin 
reference  axes  of  a  SDFI  gyroscope. 

The  existing  equation  of  motion  is  based  on  a  simplified  dynamic 
model  which  accounts  for  damping,  compliance,  and  inertia.    As  no 
closed  form  solution  to  the  equation  is  known,  a  series  form  was  used 
which  gives  a  very  clear  representation  of  the  gyroscopic  drift  terms. 
The  series  solution  appears  easy  to  handle  and  analyze  in  view  of  the 
theory  of  superposition.  The  terms  seem  reasonable  in  magnitude  and 
should  be  easily  verified  by  the  Massachusetts  Institute  of  Technology 
angular  vibrator. 
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CHAPTER  2 

THE  FLOATED  SINGLE  -DEGREE  -OF-FREEDOM 
INTEGRATING  GYROSCOPE 

2. 1      Design  Features 

The  types  and  locations  of  components  within  the  ideal  floated 
single -degree -of -freedom  integrating  gyroscope  unit  is  pictorially 
shown  in  Fig.  2-1.    These  components  are  the  gyroscopic  element, 
gimbal,  viscous  damper,  microsyn  torque  generator,  microsyn  signal 
generator,  gimbal  bearings,  and  the  case. 

The  gyroscopic  element  and  gimbal  are  hermetically  sealed  in 
a  cylindrical  shell  which  is  separated  from  the  outer  case  by  a  high 
density  fluid.    The  case  is  also  hermetically  sealed.    The  two  shells 
form  a  thin-walled  cylinder  that  performs  two  vital  functions  when 
completely  filled  with  a  high  density  fluid: 

1.  it  provides  a  viscous  damping  torque  that  opposes  the 
precessional  torque  generated  by  the  gyro  element, 

2.  it  provides  a  buoyant  force  that  neutralizes  the  weight 
of  the  gyro  element  and  its  supporting  structure. 

The  clearance  between  the  gimbal  float  and  the  case  plus  the  fluid 
viscosity  determines  the  damping  of  the  viscous  damper;  the  fluid 
density  and  the  amount  of  displaced  fluid  determine  the  buoyant  force. 


Power  is  conducted  to  the  hysteresis  type  alternating-current 
motor  by  flex  leads.    The  entire  float  assembly,  or  torque  summing 
member,  is  initially  aligned  with  the  case  by  jewel  bearings.    Final 
alignment  is  by  magnetic  suspension.    The  suspension  (elastic  re- 
straint to  float  radial  motions)  is  provided  by  the  microsyns  in  con- 
junction with  their  primary  functions  of  angle  measurement  and  torque 
generation.    The  magnetic  forces  necessary  for  suspension  originate 
with  the  change  in  excitation  currents  of  properly  tuned  microsyn 
coils  as  a  result  of  the  radial  displacement  of  the  microsyn  rotor  from 
its  centered  position.    Figure  2-2  is  a  sectioned  view  of  a  typical 
floated  single -degree -of -freedom  integrating  gyroscope.    A  more  de- 
tailed discussion  of  the  single -degree-of -freedom  integrating  gyroscope 
may  be  found  in  Refs.   1  and  2. 

2. 2     Operation 

An  ideal  single -degree -of -freedom  integrating  gyroscope  receives 
as  inputs,  the  angular  velocity  of  its  case  with  respect  to  inertial  space 
about  the  input  axis  and  electrical  current  input  signals  to  the  torque  gen- 
erator.   The  unit  generates  an  output  voltage  proportional  to  the  inputs. 
In  open  loop  operation,  the  gyro  output  signal  would  represent  the  angu- 
lar deviation  of  the  gyro  case  from  a  reference  orientation.    The  refer- 
ence is  established  from  initial  conditions  and  by  command  signals  to 
the  torque  generator. 

In  operation,  the  gyro  unit  is  housed  in  a  thermally  controlled 
environment  that  maintains  the  temperature  to  within  a  few  degrees 
Fahrenheit.     This  control  is  important  because  temperature  affects  the 
density  and  viscosity  of  the  fluid. 
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2.  THE  SYMBOL    AftrefHcrrpd)]   REPRESENTS  THE  ANGLE  A  MEASURED  FROM  THE 


EFERENCE  DIRECTION. 

COMPARED  DIRECTION 

ANGLE 

^(refMatpd)] 


DEFINITIONS: 

CASE,  (ca)  or  (e) 

THE  STRUCTURE  THAT  GIVES'SUPPORT  FOR  THE 
INTERNAL  WORKING  PARTS  OF  THE  GYRO  UNIT, 
ENCLOSES  THESE  PARTS,  AND  CARRIES  PROVI- 
SIONS FOR  EXTERNAL  CONNECTIONS  OF  ALL  KINDS. 

TORQUE  GENERATOR  AND  MAGNETIC  SUSPENSION 

UNIT,  (tg)  +  (ms) 

COMPONENT  FOR  (a)  RECEIVING  INPUT  SIGNALS  AND 
PRODUCING  A  CORRESPONDING  OUTPUT  TORQUE  AP- 
PLIED TO  THE  GIMBAL  ABOUT  THE  OUTPUT  AXIS  AND 
(b)  ALSO  PROVIDING  MAGNETIC  FORCE    ACTING  TO 
CENTRALIZE  THE  PIVOT  SHAFT  IN  THE  BEARING. 

DAMPER,  (dmp)  or  (d) 

SUBSYSTEM  RECEIVING  THE  ANGULAR  VELOCITY  OF 
THE  GIMBAL  WITH  RESPECT  TO  THE  CASE  AS  ITS  INPUT 
AND  PRODUCING  AS  THE  OUTPUT  A  RETARDING  TORQUE 
ACTING  ON  THE  GIMBAL  ABOUT  THE  OUTPUT  AXIS  WITH 
A  MAGNITUDE  PROPORTIONAL  TO  THE  MAGNITUDE  OF 
ITS  ANGULAR  VELOCITY  INPUT. 


GYRO  UNIT,  (gu) 

THE  ENTITY  MADE  UP  OF  THE  COMPONENTS  REPRESENTED 
IN  THIS  DIAGRAM  AND  ALL  THE  ADDITIONAL  PARTS  NECES- 
SARY FOR  A  SINGLE  PACKAGE  TO  CARRY  OUT  THE  FUNC- 
TIONS OF  A  GYRO  UNIT. 

SIGNAL  GENERATOR  AND  MAGNETIC  SUSPENSION  UNIT, 

(sg)  +(ms) 

COMPONENT  FOR  (a)  RECEIVING  THE  ANGLE  OF  THE  SPIN 
AXIS  WITH  RESPECT  TO  THE  CASE  AS  ITS  INPUT  AND  PRL>- 
DUCING  A  CORRESPONDING  SIGNAL  THAT  SERVES  AS  THE 
OUTPUT   SIGNAL    FROM  THE  GYRO  UNIT  AND  (b)  ALSO  PRO- 


GIMBAL,  (gim) 

STRUCTURE  CARRYING  THE  BEARINGS  FOR  THE  SPINNING 
ROTOR  OF  THE  GYROSCOPIC  ELEMENT,  ROTORS  FOR  THE 
TORQUE  GENERATOR  AND  SIGNAL  GENERATOR,  PART  OF 
THE  DAMPER,  FLOAT  SEALS  AND  STRUCTURE,  BALANCE 
ADJUSTMENTS,  STOPS,  PIVOTS,  ETC. 
RIGID  ARM  UNBALANCE  COMPENSATOR,(rouc) 

SCREW  ADJUSTMENTS  FOR  SHIFTING  THE  CENTER  OF  MASS 
OF  THE  GIMBAL  FLOAT  TO  A  POSITION  SUBSTANTIALLY  ON 
THE  OUTPUT  AXIS. 


*  A  DISCUSSION  OF  GENERALIZED  CONVENTIONS  FOR  SELF-DEFINING  SYMBOLS  OF  WHICH  Af(rof).(c 
DRAPER,  McKAY  AND  LEES  IN  INSTRUMENT  ENGINEERING  (10),  VOL.  I. 


pd)l 


IS  AN  EXAMPLE  IS  GIVEN  BY 


Fig.  2-3   Line  schematic  diagram  fer  the  single-degree-of-freed»m 
integrating  gyr»  unit. 
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The  torque  and  signal  generators  are  microsyn  units.  The  rotors 
of  these  units  are  attached  to  the  gimbal  float  while  the  stators  are  fixed 
to  the  case. 

A  line  schematic  diagram  is  shown  in  Fig.  2-3  and  a  functional 
diagram  in  Fig.   2-4.    If  the  gyroscopic  element  output  is  the  only  torque 
applied  to  the  torque -summing  member  of  the  gyro  unit,  it  is  balanced 
by  the  viscous  shear  torque.    The  viscous  torque  is  proportional  to  the 
angular  velocity  of  the  floated  gimbal  assembly  relative  to  the  case. 
This  torque  balance  causes  the  output  angular  velocity  to  be  proportional 
to  the  input  angular  velocity.    Thus,  the  angle  of  the  floated  gimbal  with 
respect  to  the  case  is  proportional  to  the  angular  displacement  of  the 
case  about  the  input  axis,  measured  with  respect  to  inertial  space.    In 
other  words ,  the  angle  through  which  the  float  turns  with  respect  to  the 
case  represents  the  time  integral  of  the  gyro  input  angular  velocity. 

The  torque  generator  permits  introduction  of  command  orienta- 
tion signals.    The  torque  generator  applies  a  torque  to  the  float,  propor- 
tional to  the  current  supplied  the  torque  generator.    If  there  is  no 
gyroscopic  torque,  the  output  of  the  torque  generator  causes  the  float 
to  rotate  (against  the  viscous  shear  torque  of  the  damper)  at  a  rate 
proportional  to  the  current  input  to  the  torque  generator.    The  gyro 
output  is  therefore  equivalent  to  an  angular  displacement  about  the  gyro 
input  axis.     The  output  axis  angular  displacement  is  measured  by  the 
microsyn  signal  generator.     The  output  voltage  of  the  signal  generator 
is  utilized  for  indication  or  control. 

For  the  gyroscope  to  be  a  linear  instrument,  the  angle  through 
which  the  float  turns  with  respect  to  the  case  must  be  kept  small.     For 
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the  integrating  gyro  this  implies  closed  loop  operation.    There  are 
two  types  of  closed  (servo)  loops  presently  in  common  use: 

1.  a  rate  feedback  loop  in  which  the  signal  generator  output 
voltage  is  amplified  and  applied  to  the  gyro  torque  gener- 
ator to  act  as  an  elastic  restraint, 

2.  a  stabilization  servo  loop  in  which  the  signal  generator 
output  voltage  is  amplified  and  applied  to  a  servo  motor 
that  supplies  an  angular  rate  about  the  gyro  input  axis. 
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CHAPTER  3 

INACCURACIES  DUE  TO  SIMULTANEOUS 
VIBRATIONS  ABOUT  THE  CASE  INPUT  AND  SPIN  AXES 

3.  1        General 

This  chapter  deals  with  the  response  of  the  single -degree -of  - 
freedom  integrating  gyroscope  when  simultaneous  angular  vibrations 
are  applied  about  the  case  input  and  spin  axes,    The  resulting  error 
torques  are  determined  and  their  physical  origin  analyzed. 

It  is  assumed  that  when  an  angular  vibration  is  present  in  the 
plane  of  the  gyro  case  input  and  spin  axes,  components  exist  about  both 
axes.    The  torque  applied  externally  to  the  gyroscope  is  developed  using 
basic  gyro  theory. 

3.  2         Coordinate  System  and  Symbols 

3.  2. 1    Coordinate  System 

The  coordinate  systems  and  the  relationship  of  the  angular  ve- 
locity components  are  shown  in  Fig.  3-1.    The  spin  reference  axis 
(SRA),  input  reference  axis  (IRA)  and  output  reference  axis  (ORA)  are 
orthogonal  axes  fixed  to  the  gyro  case.     The  gyroscopic  element  is 
considered  to  be  rigidly  fixed  to  the  float  so  that  the  float  input  axis  (IA) , 
spin  axis  (SA)  and  output  axis  (OA)  are  also  orthogonal. 
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OA 


(IC)SA 


W(IF)SA    "   W(IC)SA  cos  A(gim)   +   W(IC)IA  Sin  A(gim) 


W(IF)IA    "   W(IC)IA  cos  A(gim)   "   W(IC)SA  sln  A(gim) 


Fig.  3-1.      Gyroscope  Coordinate  System  and  Angular  Velocity 

Components 
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3.  2.  2    Definition  of  Symbols 

The  following  symbols  and  definitions  will  be  used  in  the  analysis. 


A(C-F)OA"  A(gim) 


IA 

IRA 

ORA  =  OA 

SA 
SRA 


W, 


w 


SP 


(IC)IA 


w 


(IF)IA 


W 


(IC)SA 


W 


(IF)SA 


LIA 


lSA 


Angle  of  float  with  respect  to  case 

Float  input  axis 

Case  input  axis 

Output  axis  of  both  the  case  and  float 

Spin  axis 

Spin  reference  axis  (coincident  with  SA 

when  A(gim)  "  0) 

Angular  velocity  of  gyro  wheel  about 
spin  axis 

Angular  velocity  of  case  with  respect  to 
inertial  space  about  case  input  axis 

Angular  velocity  of  float  with  respect  to 
inertial  space  about  float  input  axis 

Angular  velocity  of  case  with  respect  to 
inertial  space  about  case  spin  axis 

Angular  velocity  of  float  with  respect  to 
inertial  space  about  float  spin  axis 

Moment  of  inertia  of  wheel  and  float 
about  float  input  axis 

Moment  of  inertia  of  wheel  and  float 
about  float  spin  axis 
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Lw 


H 


;NS 


lSP 


Moment  of  inertia  of  gyro  wheel  about 
float  spin  axis 

Total  angular  momentum 

Nonspin  angular  momentum 

Angular  momentum  of  spinning 
wheel 


M 


(app) 


dH 


dt 


dH 

dt 


M 


(out) 


i,  J,  k 


V  V  nzn 

o     Jo       o 


Torque  applied  to  the  gyro  unit 
input  axis 

Rate  of  change  of  angular  momentum 
with  respect  to  inertial  space 

Rate  of  change  of  angular  momentum 
in  float  coordinates 

Torque  output  about  gyro  output  axis 

Unit  vectors  along  orthogonal  axes 
forming  right-hand  set 

Linear  acceleration 

Linear  acceleration  component  along 
platform  axis  having  i  as  unit  vector 

Linear  acceleration  component  along 
platform  axis  having  j  as  unit  vector 

Linear  acceleration  component  along 
platform  axis  having  k  as  unit  vector 

Linear  acceleration  components  along 
non-rotating  axes  having  i,  j,  k  as  unit 
vectors 


^,  0,  0 


Euler  angles 
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3.  3        Gyro  Response  to  Angular  Vibrations 

3.  3. 1    Applied  Torque  due  to  Angular  Vibrations 

In  vector  notation  the  torque  applied'  '  to  the  gyro  is 

M,ann,  =   &-       -     ^1        +   W,m  x  if  (3-1) 


dH 

dt 

I 

dH 
dt 

(app)  "  |S"  j  ■    T\  F  «F> 

By  Newton's  Law  of  Action  and  Reaction,  the  gyro  unit  output  torque  is 

dH 


M(out)=-M(app)  =  HXW(IF) 


dt 


(3-2) 


Using  the  relationships  shown  in  Fig.  3-1,  the  following  vectors  may  be 
written: 

WflF)     =  Tw(IF)SA  +  ^"W(IF)IA 

"   =  FNS  +  «SP  "  r(ISAW(IF)SA  +  ht^SP*  +T(LtAWaF)IA> 

W(IF)  X  H  =  k  [{IIA  -  ISA)  W(IF)SA  W(IF)IA  +  HSP  W(IF)IA  1      (3-3) 

The    relationships  from  Fig.  3-1  are  then  substituted  into 
Eq.  (3-3)  to  obtain  Eq.  (3-4). 

^(app)  "  ¥C<IIA-W(cos2A(gim)-slll2A(gim)^Wac>SAWffC)Li 

+    <XIA  -  W  (w2(IC)IA  "  w2aC)SA)sin  A(gim)  cos  A(gim) 

"   HSP  <W(IC)IA  cosA(gim)  "  W(IC)SA  slnA(gim)>  ]     (3"4) 
For  a  practical  single -degree -of  -freedom  gyroscope  the  gimbal  angle, 

(A) 

Ar  i    \  usually  does  not  exceed  one  degree. v  '  Therefore,  the  equation 
may  be  simplified  by  making  use  of  small  angle  approximations  and 
neglecting  second  order  terms  in  A/  .    *.    The  torque  equation  is 
reduced  to: 
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M(app)  =  k  [(IIA  "  W  W(IC)SA  W(IC)IA  +  ^IA  "  W  A(gim)    x 

2  2 

(W(IC)IA  -W(IC)SA)  "  HSP(W(IC)IA-A(gim)WaC)SA)] 

(3-5) 

For  a  single -degree -of -freedom  gyroscope  only  torques  about 

the  output  axis  are  considered.    The  rate  of  change  of  angular  momentum 

in  float  coordinates  does  not  contribute  an  output  axis  torque.    Equation 

(3-5)  represents  the  applied  torque  of  the  gyro  and,  since  it  is  all  along 

a  single  axis,  vector  notation  will  not  be  necessary. 

M(app)  "  [<XIA  "  W  W(IC)SA  W(IC)IA  +  <4a  "  W  A(gim)(W(IC)IA  "  W(IC)S^ 

"    HSP  <W(IC)IA  "  A(gim)  <W(IC)SA>  1  (3"5a> 

The  terms  of  Eq.   (3-5)  represent  error  or  undesired  torques. 
A  single-degree-of-freedom  gyroscope,  operating  in  a  closed  loop, 
would  have  an  angular  displacement  from  the  null  or  zero  position  that 
would  be  in  error  due  to  these  undesired  torques.    An  examination  of 
these  error  torques  will  now  be  made. 

3.  3.  2    Dynamic  Torques 

The  first  torque  component  to  be  considered  will  be 

»IA-WW(IC)SAWaC)IA  <3-6> 

It  is  evident  that  the  nature  of  the  error  will  depend  upon  the  variation 

of  the  two  simultaneous  angular  velocities  and  the  float  inertias  about 

(5) 
the  spin  and  input  axes.    This  error  is  known  as  anisoinertiav  '  or 

dynamic  unbalance.    It  could  be  eliminated  by  building  gyro  units 

having  equal  inertias  about  both  axes. 
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Next  consider  the  second  term  of  Eq.  (3-5) 

«IA  "  W  A(gim)  (W(IC)IA  "  W(I(W  <3"7) 

In  the  above  expression  the  nature  of  the  error  will  be  determined  by 
the  difference  of  the  inertias  as  well  as  the  difference  of  squared  angu- 
lar velocities.    Variation  of  gimbal  angle  must  also  be  considered  since 
it  is  not  a  constant. 

Expressions  (3-6)  and  (3-7)  are  dynamic  torques^  '  and  may  be 
interpreted  as  representing  four  angular  momentums  and  their  associa- 
ted precessional  angular  velocities.    Consider  that  there  are  four  angular 
velocity  components  and  four  corresponding  angular  momentums. 

Momentum  Components  Angular  Velocity  Components 


%AWaC)SA 

W(IC)SA 

-A(gim)IIAW(IC)SA 

~A(gim)W(IC)SA 

IIAW(IC)IA 

W(IC)IA 

A(gim)ISAW(IC)IA 

A(gim)W(IC)IA 

The  total  "dynamical  torque'"  '  when  each  forced  precession  is 
considered  is: 

MD      =    "  *SA  W(IC)SA(W(IC)IA  "  A(gim)  W(IC)SA^ 

-^A^gim)  W(IC)SA(W(IC)SA  +  A(gim)  W(IC)IA) 

+  ^A^lOlA^dOSA-"  A(gim)W(IC)IA) 

~  %  A  A(gim)  W(IC)IA  (W(IC)IA  "  A(gim)  W(IC)SA^ 
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The  part  in  parentheses  for  each  of  the  above  terms  is  the  forced  pre- 

cessional  rate  while  the  part  before  the  parentheses  is  the  angular 

momentum.    For  example,  the  angular  momentum  component,  IsAWfIC)SA 

precesses  about  the  input  axis  with  the  angular  velocity, 

W(IC)IA  ~  A(gim)W(IC)SA°    When  tne  equation  is  expanded  and  second 

order  terms  of  A/  .    %  neglected,  the  expression  for  dynamic  torque  is 

obtained. 

MD  =  tflA-W  W(IC)IAW(IC)SA  +  "iA^Sa' A(gim)(W(IOIA- W(2IC)SA^ 

(3-8) 
3.  3.  3    Gyroscopic  Torque 

The  applied  gyroscopic  torque  on  the  gyro  is 

Mge    =    -HSP(WaC)IA-A(gim)W(IC)SA)    =-HSPW(lF)IA 

(3-9) 

Variation  of  this  torque  will,  as  in  previous  considerations,  be  de- 
pendent on  the  angular  velocities  and  the  gimbal  angle.    The  component, 

HSP  A(gim)  W(IC)SA  '  is  of  sPecial  significance  and  is  called  Coning 
Error.  ^' 

3. 4        Coning 

3.  4. 1    General  Discussion  of  Coning 

The  coning  error  defined  above  is  more  usually  considered  a 
cross  coupling  error.    It  is  a  result  of  oscillations  in  the  plane  of  the 
case  input  and  spin  axes.    Comparing  the  form  of  the  coning  expression 
from  the  next  section  with  that  of  cross  coupling  reveals  a  close  similar- 
ity.   Because  of  this  similarity,  it  is  possible  to  consider  cross  coupling 
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as  another  cause  of  the  coning  error.    This  form  of  coning  error  may 
be  reduced  by  restraining  rotations  about  the  output  axis.    A  tight  servo 
loop  would  keep  the  error  negligible. 

If  the  gyro  input  axis  has  a  motion  such  that  it  returns  periodically 
to  its  original  position,  then  the  input  axis  will  trace  out  a  path  in  space. 
This  motion  is  known  as  the  coning  motion  and  is  present  even  in  a  theo- 
retically perfect  instrument. 

The  coning  error  is  a  geometric  effect  that  allows  the  gyro  to 

(7) 
respond  to  an  input  rate  produced  by  the  coning  motion.         Almost  any 

motion  of  the  case    creates  coning.    Consider  oscillations  in  the  plane 

of  the  case  input  and  spin  axes.    The  component  along  the  input  axis 

will  cause  the  gyro  to  precess  about  the  output  axis.    When  the  oscillating 

motion  about  the  spin  reference  axis  is  combined  with  this  precession,  the 

result  is  a  coning  motion  of  the  float  input  axis.    As  another  example,  if 

there  is  oscillation  of  the  case  in  the  spin  axis  -  output  axis  plane,  coning 

would  again  result  because  the  damping  fluid  will  drag  the  float  around. 

When  single -degree -of -freedom  gyros  are  used  for  stabilization, 
a  coning  error  exists  that  cannot  be  eliminated. 

3.  4.  2    Platform  Coning  Effect  ^8^ 

Consider  two  coordinate  systems: 

i  ,  j   ,  k  platform  coordinate  system 

y    y     y 

i  ,  j      k  coordinate  system  defining  the 

orientation  the  platform  servos 
are  trying  to  maintain 
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ko>kl 


* 


jp  j2 


Fig.   3-2.      Euler  Angle  Relations  Between  Platform  and  Desired 

Orientation  Coordinates 
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By  use  of  the  Euler  angle    rotations5taken  in  the  order  shown  in 
Figo   3-23  the  coordinate  transformations^   using  matrix  notation,  are  as 
follows: 


and 


and 


cos  \p- 

-sin  yp 

0 

sin  \p 

cos  \p 

0 

0 

0 

1 

w10    =   ^ko   =  ^kl 


Tl 

cos  e 

h 

= 

0 

ki 

-sin  9 

sin  9 

T2 

0 

T2 

COS0 

*2 



'21 


*h    =    ^2 


]2 
^2 


0 

COS  0 

sin  0 


0 

i 
P 

sin  0 

rF 

COS  0 

k 
P 

— 

to). 


•  — 
0  i 


(3-10) 
(3-11) 


(3-12) 
(3-13) 


'p2    -    ^  x2         *    p 

It  is  possible  to  combine  the  three  matrices  and  simplify  them  by 
making  the  standard  small  angle  approximations  : 


(3-14) 
(3-15) 


cos  \p    -sin  ip     0 

sin  >p       cos  \p    0 

0  0         1 


cos  9     0      sin  9 

0  1        0 

■  sin  9      0      cos  9 


1  0  0 

0      cos  0     -  sin  0 
0      sin  0         cos  0 


ft 
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1 


1 


(3-16) 


The  angular  velocity  of  the  platform  coordinate  system  with 
respect  to  the  orientation  coordinate  system  may  be  defined  as: 


w 


po 


pip  +  qjp  +  rkp 


(3-17) 


From  Eqs.   (3-10),  (3-11),  and  (3-12)  we  also  have 


po 


"    V+a,21  +  a,10 


=    0i    +  0  j2  +  tf/kj  (3-18) 

In  order  to  express  Eq.   (3-18)  in  platform  coordinates,  substitute 
relation  for  k.,   from  Eq.  (3-12)  and  then  make  use  of  Eq.   (3-14)  as 


follows: 


po 


$  i    +  5  j9  +  ^(nsin  0  i9  +  cos  0  k9) 


J  i    +   [-  if  sin  0      5       if     cos  0  ] 


^2 


=  $  i    +  [-if  sin  0      0     i? cos  0] 


1        0  0 

0    cos  0    -sin  0 
0    sin  0     cos  0 


=  ($  -  if  sin  0  )  i    +  (0  cos  0  +  if  cos  0  sin  0)  j 

+  (-0  sin  0  +  if  cos  0  cos  0)  k 

Assuming  small  angles,  a  comparison  with  Eq.   (3-17)  gives 


23 


p   =   0  -e$ 

q    =     0    +  0  >p 

r    =     '  -0^  (3-19) 

The  rates  of  change  of  the  Eulerian  angles  are 

0  =  p+6p  -  p+(9(r  +  0  0)=rp+r# 
0  =  q-01//  =  q-0(r  +  0  0srq-r0 
'/>    =    r+0#     =    r  +  0(q-0^)   *  r  +  0q  (3-20) 

Integrating  Eq.   (3-20)  to  obtain  the  expressions  for  the  Euler  angles 

we  obtain: 

0(t)     =     /  pdt  +    /  r9  dt 
o  o 

t  t 

6  (t)     =    /  q  dt  -     /  r0    dt 
o  o 

t  t 

\{/{t)     =     f  r  dt  +     f  q0   dt  (3-21) 

o  o 

In  the  above  equations,  0 ,   9 ,  and  *//  are  angles  in  space  and  when  all 
are  zero,  the  platform  coordinate  system  coincides  with  the  orientation 
coordinate  system  by  definition.    Having  the  space  angles  zero  is  no 
guarantee  that  the  integrals  of  p,  q  and  r  are  zero.    Depending  upon  the 
past  history  of  motion,  p,  q  and  r  may  or  may  not  be  zero.     For  ex- 
ample, having   0(t)  =  0   does  not  ensure  that    /    0rdt   is  zero.    A 

o 
typical  platform  stabilization  system  can  maintain  the  integrals  of  p, 

q  and  r  at  zero.    Therefore,  the  error  of  the  platform  coordinate  sys- 
tem is  given  from  Eq.  (3-21)  by 
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t 

4  0  (t)      =       /    e  r  dt 
o 

t 
A  0(t)      =    -  f  0r  dt 
o 

t 

A^(t)       =         /0qdt  (3-22) 

0 

If  the  vibratory  motions  are  small,  the  integrals  above  become 
significant  only  if  there  is  a  rectification  of  the  motions.    Let 
9    -    6      sin  art 

r    =    r      sin  (art  -  y ) 
then 

A  0  (t)    -    Or      f  sin  cot  (sin  art  cos  y  -  cos  art  sin  y  )  dt 
o 

r1         2 
=  0     r        /  (sin    art  cos  y-sin  a?t  cos  art  sin  y)  dt 

o 

e  „  r    r 


o    o 


—  sin  2art  |cos  y  +  —  (1  +  cos  2art)  sin  y 

4  /  4 


art 
2 

(3-23) 


An  examination  of  Eq„  (3-23)  reveals  one  linear  term.    Therefore,  the 
angle  drifts  from  null  with  a  constant  drift  rate. 

9„rn  cos  y 

a  0(0   =  _^_2 

2 

For  the  system  just  analyzed,  there  will  exist  a  coning  error  that 
grows  with  time.    One  possible  method  of  reducing  the  error  would  be  by 
use  of  a  coning  computer, v  ;    The  computer  would  solve  Eqs.   (3-20)  and 
compensation  for  the  coning  drift  could  then  be  built  into  the  system. 
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3.  5        Sculling  Effect 

Sculling  produces  an  acceleration  error  in  accelerometers 
mounted  in  gyro  stabilized  platforms  undergoing  combined  linear  and 
angular  motion.    This   effect  will  be  considered  because  what  is  really 
important  is  the  effect  of  this  error  on  position  indication. 

The  same  coordinate  systems  and  coordinate  transformation 
matrix  previously  used  in  describing  platform  coning    apply.    Define 
the  total  acceleration  vector  in  both  coordinate  systems  as 

"   =   nx  fp  +  ny  ip  +  nz  *p  (3"24) 

°   "   nxnro  +  nyJo+V/o  <3-25> 

o  Jo  o 

If  there  is  no  acceleration  error,  then 


nx 

— 

\ 

- 

0 

ny 

- 

X 

= 

0 

nz 

- 

X 

= 

0 

(3-26) 

A  check  on  the  above  relations  will  now  be  made.     From  transformation 
(3-16)  we  obtain 

n   =   nxJIp  -  ^Tp+  9  ip]   +    n^  [  ^Tp +Tp  -  0  kp] 

+  nz    [-  9  Tp+0Tp  +  "kpI      (3-26) 
Comparison  of  Eqs.   (3-24)  and  (3-26)  produces  the  relations 
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■*- V*"V\ 


n     =    ~^n      +  n      +  0n 
*  o       yo  o 


n*  =  8\-*X  +  X 


and  by  rearranging 


nx  ~  nx     =    ^nv    "  e  nz 
x        xo  yo  zo 

n    -  n       =    0  n      -  \Z/n 
y        yo  zo  xo 

nz  -  nz      =    0  nx    "  *  nv  (3_27) 

o  o  yo 

Equations  (3-27),  when  compared  to  Eqs.   (3-26),  show  that  errors  be- 
tween the  accelerations  actually  measured  and  the  accelerations  we 
desire  to  measure  exist. 

Assume  a  component  of  the  vibration  of  the  platform  is  at  a 

certain  frequency,   oj.    Then  consider  any  term  of  Eq.   (3-27),  say 

0  n      .    Let  the  linear  and  angular  vibration  vary  sinusoidally  at 
o 


frequency  w. 

If  0    =  |  0  I  sin  cot 


nz    =  |nz  I  sin(wt  -  y) 


then 


0n        =    1 0  I  In.     I   sin  cot  (sin  cut  -  y  ) 
zo        '    "  zo 


nzm 

— _  [  ( 1  -  cos  2o)t)  cos  y  -  sin  2cut  sin  y  ] 
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which  gives  a  steady  state  drift  of 

|0  nr 


zol 


COS  y 


2  (3-28) 

As  a  consequence  of  the  analysis,  it  is  seen  that  small  angular 
motions  of  the  platform  can  result  in  rectified  acceleration  drifts  whose 
magnitudes  depend  on  the  phase  relationships  between  the  linear  and 
angular  motions.    The  name  "sculling"  v  '  is  given  to  this  effect  because 
the  motion  causing  the  drift  is  identical  to  that  produced  in  propelling  a 
boat  with  an  oar  at  the  stern. 
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CHAPTER  4 
MAGNITUDE  OF  THE  ANGULAR  VIBRATION  INACCURACIES 


4. 1        General 

The  errors  producing  forced  precession  of  a  single-degree-of- 
freedom  gyroscope  have  been  analyzed  in  the  previous  chapter.    This 
chapter  deals  with  a  more  practical  consideration;  the  magnitudes  of 
the  various  error  terms.    When  practical,oscillations  capable  of  being 
produced  by  the  Angular  Vibrator ,  built  by  the  Instrumentation  Labora- 
tory of  the  Massachusetts  Institute  of  Technology,  will  be  assumed. 

The  following  are  assumed  as  typical  values  for  a  single -degree  - 
of -freedom  integrating  gyroscope: 

2 


A I    =   Iia-IsA    =    ~  65°  gm  cm' 

H    -   HQT5    =    2  x  106  gmcm 
SP  sec 


=    1 


Cd 


A/  •    \    <    •  01745  radians 

A  comparison  of  the  various  inaccuracies  previously  considered 
will  now  be  made. 
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4.  2        Magnitude  Determination 

For  convenience,  the  inaccuracies  developed  in  the  previous 
chapter  are  again  listed. 

MD  =  (IIA-ISA>W(IC)IAW(IC)SA+  ^W  A(gim)(W(IC)IA-  W(IC)SA> 

(3-8) 

Mge  =    "HSP  (W(IC)IA  "  A(gim)  W(IC)SA>  <3-9> 

t 
A  0  (t)    =  /   6  r  dt 
o 
t 
A0(t)  =  -/0rdt 
o 
t 
A<Mt)  =   /     0qdt  (3-22) 

o 

nz  "  nz„   ■    9  nx„  -  *  V  (3-27) 

o  0  Jo  v  ' 

Equations  (3-8)  and  (3-9)  apply  to  a  single  gyroscope  instrument  while 
Eqs.  (3-22)  and  (3-27)  apply  to  a  platform  stabilized  system.    The  mag- 
nitudes of  the  various  terms  will  be  determined  for  the  case  of  harmonic 
angular  velocities. 

4.  2. 1    Magnitude  of  the  Dynamic  Torque 

Let 

W(IC)IA   =   !wIA|sin(wit+0i> 

W(IC)SA    =|WSA!Sin(w2t+^2>  <4"» 
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The  anisoinertia  component  of  Eq.   (3-8),  when  using  Eq.   (4-1) 
may  be  written 

W, 


Iai 


WIA 


gA       COS  [(U1-U)2)t+  (01-02)1 

COS  [  (Wj  +  u>2)t  +    0i+02l  (4~2) 


wl  =  w2  =  w 


If 

then  the  above  reduces  to: 

—  |WIa|     WSa|     Icos(  0'l  "  ^  "  cos(2wt+  0i  +  02>1 

(4-3) 

Equation  (4-3)  consists  of  a  constant  term,  and  a  cosine  term  that 
has  a  zero  average  value.    Since  (0  .  -  0  «)  represents  the  phase  dif- 
ference between  the  two  velocities  it  is  evident  that  only  the  in  phase 
components  of  the  motions  create  an  error.    If  the  two  angular  veloci- 
ties are  in  quadrature,  there  will  be  no  resultant  error. 

If  b)«  \  o>2  5  then  the  torque  will  oscillate  at  a  beat  frequency 

(cjj  +  aO  ,    (wj  -  co2)  and  have  an  average  value  of  zero. 

The  second  component  of  Eq.   (3-8)  for  u)  =  o><  =  a)2  may  be 
written  as: 

f-  A(gim)  [   IWIa!2  "IWSa|2  +  cos  2^(|WSA|  2  cos  202 

+  |WIA|  2  cos  20  x)  -  sin2wt(|WSAl  2  sin  20  2 

+  |WIAl2sin201)]  (4-4) 
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Before  the  magnitude  can  be  determined,  the  variation  of  gimbal 
angle  with  time  must  be  known.     For  the  present,  an  ideal  integrating 
gyro  will  be  assumed  so  that 

A(gim)    =   W(IC)IA 
and 

A^  -  I  wnniA  dt 


(gim)    '     J    "(IQIA 

-IwiaI 


(J) 

Substituting  (4-5)  into  (4-4)  we  get 


cos  (wt  +  0X)  ^45j 


~—   I W J    sin  ( cot  +  0  x)  F|WIA|  2  -  |ws J 
2a;  L_ 


2i 

+  cos  2a>t(|wSAl  2  cos  20  2  +  IWjJ  2  cos  20  j) 

-  sin  2o>t(|wSA|  2  sin  20  2  +   (Wj^l  2  sin  20  Jj  (4-6) 

AH  the  terms  in  Eq.   (4-6)  have  zero  as  an  average  value;  therefore, 
the  entire  expression  may  be  neglected.    The  same  applies  for  the 
case  of  o>]^  u>2* 
4.  2.  2    Magnitude  of  the  Gyroscopic  Torque 

We  next  turn  our  attention  to  Eq.   (3-9).    After  making  substi- 
tutions for  the  motions  we  have 

Mge  =  "  HSP[IWI^    sin  (<"!*+  #1>  -  A(gim)(wsi  sin("2t+  #2> 

(4-7) 

The  first  term  in  the  bracket  will  have  a  zero  average  value.    Our 
attention  will  now  focus  on  the  cross  coupling  coning  term: 
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A(gim)  HSP  IWSa'    si«(^t  +  *2)  (4"8> 

Again  assuming  a  perfect  gyro  so  that  relation  (4-5)  applies  we  obtain 


Hsp  |wiA|   |wgAl   cos(Wlt+  0X)  sin(o>2t+  02)  (4-9) 


If 

we  can  write: 


»1  -  «2  "  w 


2w 


—  HSP  iWIA|    rsJ   tsin^2  -  *1>  +  sin(2wt  +  02  +  0'1>  I 

(4-10) 

Thus  the  torque  depends  on  the  out  of  phase  components  of  the  veloci- 
ties.   Should  the  two  velocities  be  in  phase  there  would  be  no  torque. 
Only  velocities  in  quadrature  produce  a  coning  error.    From  this 
discussion  it  is  seen  that  the  float  drifts  from  null  at  a  constant  rate 
of 

— KJ     KJ     sin(02_0l)     rad/sec  (4"n> 

2  co 

For  the  condition  of  o>*  ^  a^,  Eq.   (4-9)  contains  no  linear  or 
average  torque  values. 

Note  that  when  coning  is  a  maximum,  the  anisoinertia  effect  is 
a  minimum.    Angular  motions  that  are  not  exactly  in  or  out  of  phase 
will  produce  both  effects.    However,  even  for  this  condition  it  is  seen 
that  the  anisoinertia  error  may  be  neglected,  since  it  is  very  small 
in  comparison  to  the  cross  coupling  coning  error.    This  is  due  to  the 
fact  that  the  coning  term  involves  Hgp  which  is  very  much  larger  than 
any  of  the  values  in  Eq.  (4-3). 

The  coning  drift  rate  is  plotted  in  Fig.  4-1  as  a  function  of 
frequency  for  angular  velocities  30,  60  and  90  degrees  out  of  phase. 
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All  of  the  previous  work  has  been  for  an  ideal  gyroscope.    Ac- 
tually, the  true  gimbal  angle  variation  should  be  obtained  from  the  gyro 
equation 

A(gim)  +  W(IC)SA  A(gim)  =  W(IC)IA  (4"12> 

If  Eq.   (4-12)  is  used  to  obtain  the  gimbal  angle  motion,  the  only  differ- 
ence from  the  previous  result  would  be  a  nonzero  average  value  of  coning 
drift.    This  can  be  easily  shown.     First  integrate  Eq.   (4-12) 

A#       ^    _    e-/W(IC)SAdt      ft  e/W(IC)gAdt 

A(gim)    "    e  J0  W(IC)IA   e  dt 

(4-13) 

Assuming  sinusoidal  variations  and  small  magnitudes  for  the  angular 

velocities,  the  exponentials  may  be  approximated  by  the  first  two  terms 

of  a  series 

±/W(IC)SAdt       -  4     rw  M 

e  *  1  ±  f  W(IC)gA  dt 

(4-14) 
Using  the  relations  (4-1)  and  the  gimbal  angle  expression  (4-13),  we 
can  determine  the  nature  of  the  coning  term 

t         ,  1 

A(gim)    *      f  IwiaI    !WSa'    — —  sin(wit+  0i)  cos(ct)2t+  02)  dt 
o  a)2 

*     —   IwJ    lwgAl  i   /[|sin(w1t  +  01-w2t-02] 
it 

+  sin[o)..t  +  a>2t  +  0  ^  +  0 2]!  dt 

if  u).  =  a?2  =  oj 

then 

..                                                               cos(0  1  +  0  0)  -  cos(2cut+01  +  09) 
Atam)^lWIAl    lWSAl      sinC^-^Jxt, l-^— 1— ± 


2w  2a> 
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It  is  seen  that  a  constant  drift  rate  appears  similar  to  expression 
(4-11)  plus  an  oscillatory  term  and  a  nonzero  average  term.    Since 
the  constant  drift  rate  term  dominates,  the  approximation  for  gimbal 
angle  motion  using  Eq.   (4-5)  is  valid. 

4.  2.  3    Magnitude  of  Platform  Coning 

From  Chapter  3  it  was  seen  that  platform  coning  effect  causes 

the  Euler  angles  to  drift  from  null  at  a  rate  equivalent  to     0(t) 

0     r 
A0(t)    =        °     °     COS   y  (4-15) 

2 

The  above  equation  was  the  result  of  integrating  one  of  the  error  equa- 
tions (3-22). 

If  it  is  assumed  that 

r    =    r     sin  (o>t  +  y )  «  \p 


xp    =  p 

0    C0S 

(OJt 

+  r) 

then 

• 

i,  =  _ 

*ow 

sin  (art  + 

y) 

ro   "    - 

^ow 

and  Eq. 

(4-15) 

is 

now 

A*(t) 

- 

^0 

e0  a, 

cos 

y 

(4-16) 
2 

The  assumption  will  now  be  made  that  the  amplitudes  of  the 
motions  are  known  but  that  the  directions  and  phase  relationships  be- 
tween the  components  of  the  motion  are  random. 

Let 

\p     =   K  cos  j3 

9      -    K  sin  p 
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Hence  K2 

A0(t)    =     -    Z     sin  2/3     cosy 
4 

Now  assume  that  /3  and  y  are  independent  and  uniformly  distributed 

from  0  to  27T,  so  that  the  mean  value  of  each  of  the  above  trigonometric 

terms  is  zero,  and  the  variance  of  each  is  1/2.    For  example 

x  2ir 

Mean  (cos  y)    =  u.  cos  y    =  E(cos  y)  =  f  cos  7   dr    =   0 

2ir       o 

2                             2            1       J"'    2 
Variance  (cos  y)  =  a    (Cos  y)  =  E  (cos  y)  = J  cos   y  d.y  =    1/2 

2tt     o 
We  can  now  write 

E[A?(t)]  I)       E[  sin2/3    cos  yf 


J*2**)        r  i.      1 

\    4   /  2       2 


2 

=    a 


and  finally  the  standard  deviation  or  RMS  value  is 

a  =   K   £  (4-17) 

8 

where 

K    =    maximum  amplitude,  radians 

a>   =   frequency,  rad/sec 

Equation  (4-17)  is  plotted  in  Fig.   4-2. 

4.  2.  4    Magnitude  of  Sculling  Effect 

The  sculling  effect  will  be  treated  by  the  method  used  in  the 
previous  section.    Consider  the  first  of  Eqs.   (3-27)  : 

nx-v  =  *"v  -ev  <4-18> 

o  Jo  o 
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Using  the  form  of  the  sculling  error  found  in  Chapter  3  we  can  write 

nx"nxo  =  7  K\J  cosri  "7  KnzJ  cos^2 


(4-19) 


If  we  let 


n  =   n  cos£ 

yo 

n      |     =   n  sin  0 
zo 

|  \J/    I     =  K  cos  a 


6     I     =Ksin   a 


then 

n„  -  n       =  —  nK  (cos  a  cos  j8   cos  y<    -  sin  a  sin  /3   cos  y«) 
x        xo        2  x  ^ 

(4-20) 

Consider  that  a,  /S,   y-,  and  y«  are  independent  of  each  other  and 
uniformly  distributed  from  0  to  2ir.    The  mean  value  of  each  trigono- 
metric function  is  zero  and  the  variance  of  each  is  1/2.    Then 

212  222  222 

E  [(nx  -nx  v    ]    =    (—   nK)      E  [  cos    a  cos    /3  cos    y^  +  sin    a  sin   )3  cos   y2] 

o  ^ 

1/2 
aK-nx)    =    inK[2(-i)  (l)(i)] 
x        xo  2  2        2        2 

=    J*         standard  deviation  or  RMS  (4-20) 

4 

where  o 

n   =   maximum  acceleration  (ft/sec  ) 

K   =   maximum  amplitude  (radians) 
Equation  (4-20)  is  plotted  in  Fig.   4-3. 
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CHAPTER  5 

THE  RESPONSE  OF  THE  SINGLE -DEGREE -OF-FREEDOM 
INTEGRATING  GYROSCOPE  TO  TWO  ANGULAR  VIBRATIONS 
SIMULTANEOUSLY  APPLIED  ABOUT  THE  SPIN  REFERENCE 
AND  OUTPUT  AXES 


5.  1     Introduction 

Gyroscopic  effects  caused  by  sinusoidal  angular  vibrations  are  studied 
in  this  chapter.     Two  vibrations  of  different  frequencies  have  been  applied  to 
two  of  the  body  fixed  axes  of  the  gyroscope  case.    The  particular  axes  to  which 
these  inputs  are  applied  are  the  output  and  the  spin  reference  axis. 

A  series  solution  has  been  used  to  approximate  the  theoretical  angle  which 
should  exist  between  the  float  and  inertial  space  about  the  output  axis.    The 
method  used  was  applied  to  a  case  of  two  equal  input  frequencies  in  ref.   10. 
Thus,  this  is  the  more  general  situation  and  a  further  study  of  the  problem  with 
varying  inputs. 

A  typical  dynamic  model  of  a  gyroscope  has  been  used  as  the  basis  for 
the  equation  of  motion.    Effects  such  as  inertia  damping,  and  compliance  have 
been  included.    The  resulting  equation  presently  has  no  closed  form  solution 
and  the  approximate  series  method  was  used. 

The  investigation  was  an  attempt  to  discover  if  and  when  gyroscopic  drift 
can  be  expected  along  with  the  associated  values  of  magnitude.  . 
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5.  2     Definition  of  Symbols 

The  following  symbols  and  definitions  have  been  used  in  the  analyses: 
(units  are  in  parentheses) 

a  =  constant  (radian) 

b  =  constant  (radian) 

c  =  constant  (radian) 

cd(OA)  =  damping  coefficient  between  float  and  case  about  the 


h    = 


output  axis  (dyne-cm-sec/radian) 
H 


cd(OA) 


k    =     K(C-F)OA 

cd(OA) 
t    =    time  (seconds) 
x    =    the  variable  representing  A/j  f\q  a 

^(I-C)NS     =    anSu^ar  rotation  of  the  case  about  a  stationary  axis 

aligned  with  the  earth's  axis  which  coincides  with  the 
spin  reference  axis  when   Aq  r\QA  =  0  (radian). 

A/j  f\qa    =    angular  rotation  of  the  float  about  the  output  axis  (radian) 

A/,  q\qa    =    angular  rotation  of  the  case  about  the  output  axis  (radian) 

A  =  constant  (radian) 

B  =  constant  (radian) 

C  =  constant  (radian) 

D  =  constant  (radian) 

E  =  constant  (radian) 

F  =  constant  (radian) 

G  =  constant  (radian) 
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H    =    angular  momentum  of  the  wheel  (dyne -cm -sec) 
Hj    =    constant  (radian) 

I(F)OA    =    inertia  of  the  float  about  the  output  axis  (gm-cm  ) 
1^   =    inertia  of  the  wheel  about  a  radial  axis  (gm-cm^ 

J     =    constant  (radian) 
K    =    constant  (radian) 

K(C-F)OA  =  elastic  restraint  about  the  output  axis  between  float  and 
case.  Rotational  stiffness  between  float  and  case  about 
the  output  axis  (dyne-cm/rad) 

L  =  constant  (radian) 

M  =  constant  (radian) 

N  =  constant  (radian) 

a)..  =  circular  frequency  of  vibrational  input  (rad/sec) 

W9  =  circular  frequency  of  vibrational  input  (rad/sec) 

r     =      '    I — -     =   gyroscope  time  constant  (seconds) 

cd(OA) 

(The  inertia  of  the  wheel  is  added  to  that  of  the  float  in  t 
since  the  wheel  is  assumed  to  be  perfectly  rigid  with 
respect  to  the  float  about  the  output  axis  in  the  analyses.) 

5.  3      The  Equation  of  Motion  of  the  Dynamic  Model 

The  simplified  equation  of  motion  of  the  dynamic  model  shown  in  Fig.   5-1 
is  given  by  : 

A/ 


K(C-F)OA 


(I-F)OA  ~  A(I-C)OAj    +    cd(OA)  [A(I-F)OA  _  A(I-C)Oa| 
(I-F)OA    =    H  I  A(I-C)NS  A(I-F)OA     (5-1) 


J(F)OA  +  JW 
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[  VC)NS  A(I-F)Oa]    is  introduced 


By  using  small  angle  theory,  the  term    H  I  A,T  «*„„  A, 

by  the  component  of  angular  velocity  about  the  input  axis  produced  by  the  motion 

about  the  space-fixed  axis.    Division  by   c^q**  and  rearrangement  yields: 


K 


(C-F)OA 
'd(OA) 


A(I-F)OA  +  A(I-F)OA   + 


I(F)OA  +  *W 
cd(OA) 


L(I-F)OA 


H 


_cd(OA)  _ 


A(I-C)NS  A(I-F)OA 


K 


(C-F)OA 


'd(OA) 


A(I-C)OA  +    A(I-C)OA 


or 


k_hA(I-C)NS      A(I-F)OA    +    A(I-F)OA   +    T  A(I-F)OA 


"   A(I-C)OA   +    kA(I-C)OA 


(5-la) 


The  stiffnesses  between  the  wheel  and  float  about  the  input  and  output  axes  are 
assumed  to  be  infinite  compared  to  the  inertia  of  the  wheel  for  the  frequency 
range  considered.     The  stiffness  and/or  damping  between  the  float  and  case 
about  the  input  axis  is  also  assumed  infinite  compared  to  the  inertia  of  the 
float  and  wheel  within  the  frequency  range  considered. 


5.  4     Method  of  Solution  of  the  Equations  of  Motion  for  Two  Different 
Sinusoidal  Vibration  Inputs  of  Small  Magnitude 


The  inputs  introduced  to  the  gyroscope  case  may  be  represented  as  follows: 

(5-2) 
(5-3) 


10 


A(I-C)NS    =   a  Sin  "l* 


A(I-C)OA    =  b  sin  CJ2t  +  C  C0S  CJ2t 

In  most  applications  of  the  single-degree-of-freedom  integrating  gyroscope, 
the  constant  h  is  designed  equal  to  unity.     For  purposes  of  simplification,  this 
constant  will  be  included  in  the  vibration  inputs  as  a  factor  of  the  constant  a. 


44 


s 


il? 

♦ 

U- 

z 

U) 

0 

<? 

LU 

_J 

& 

<£ 

£ 

o 

u 

o 

U- 

o 

ijj 

o 

£ 

a 

< 

q 

o 

z 

u 

■z 

> 

(/) 

< 

P 

> 

z 

<c 

cJ5 

eft 

H 

LLi 

U- 

X 

CO 

K 

i/) 

H 

, 

3 

0 

1 

45 


To  be  strictly  correct,  earth  rate  should  be  included  in  the  definition  of  the 
vibration  input  about  the  North-South  axis.    This  rate  is,  however,  small 
compared  to  the  input  angular  velocity. 

Substitution  of  equations  (5-2)  and  (5-3)  into  equation  (5- la)  yields  a 
linear  differential  equation  with  a  variable  coefficient  of  the  following  form: 

k  -  aw.,  cos  00 A      x+  x+  rx 

=  ba>2  cos  Wgt  -  co>2  sin  a^t  +  kb  sin  ajot  +  kc  cos  aj«t  (5-4) 

where   A(I_F)QA   =   x 

There  is  no  known  closed  form  solution  for  equation  (5-4).    However,  if  the 
constants  are  small  compared  to  unity,  a  suitable  series  solution  may  be  ob- 
tained. * 

5.4.1  Negligible  Elastic  Restraint  Between  the  Float  and  Case 
About  the  Output  Axis 

When  the  elastic  restraint  is  very  small  compared  to  the  damping  con- 
stant, i.  e.  , 

k*0 

equation  (5-4)  reduces  to: 

(-acjj  cos  u)A)x+x+  tx  =   bo>2  cos  u>2t  -  ca>2  sin  a>2t  (5-5) 

Assume  a  solution  of  the  form: 

x   =   Xj  +  x2  (5-6) 

where  x.,  is  defined  by: 

xl  +    TX1    =  ^2  cos  w2*  ~  ca,2  sin  ^2*  (5-7) 

Substitute  (5-6)  and  (5-7)  into  (5-5): 

(-awj  cosc4);.t)(x1  +  x«)  +  (x.  +  x2)  +   r{x^  +  x2)    =  bcjg  cos  c^t-  co>2  sinct>2t 

(5-8) 


♦The  series  solution  technique  was  suggested  by  Mr.  Daniel  Goldenberg 
of  the  Mathematics  Group,  MIT  Instrumentation  Laboratory. 
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or 

(-aojj  cos  (jjtjxg  +  x„  +    rx2    =   (aw.cosw.t)  x1  (5-8a) 

Now  let  x2  =  x3  +  x4  (5-9) 

where  x„  is  defined  by: 

x«  +  TXo    =    (aojj  cos  w«t)  x-  (5-10) 

Substitute  (5-9)  and  (5-10)  into  (5-8a): 

(-aa><  cos  u)A)(x^  +  xj  +  (x„  +  x.)  +  t(x"    +  x\)   =  (aw.,  cos  wA)  x.  (5-11) 

or 

(-aoj^  cos  qjA)x.  +  x.  +  tx,  =  {a.u)^cosu)A)x^  (5-lla) 

Continuation  of  this  process  results  in  a  series  solution  of  the  following  form: 

x  =  xl  +  x3  +  x5  +  ••••••  •  x2m-l  +  x2m  <5"12) 

where  x«      is  the  even  term  remaining  after  the  above  process  is  repeated 
m  time  So    x2      satisfies  the  equation: 

(-ac^cos  «1t)x2m  +  x2m  +  rx2m    =    (wlCosuA)x2ml  (5-13) 

A  Brief  Summary  for  k«*  0  is: 
Problem: 

(-aw-,  cos  (j)+  t)  x  +  x  +  tx   =   bco2  cos  <joA  -  ca>2  sina>2t  (5-5) 


let 
where: 


x  =  x1  +  x3  +  x5+.    .   .   .    x2m_1  +  x2m  (5-12) 


x<  +  tx-     =    ba)2  cosa>2t -ca>2  sinajgt  (5-7) 

L  +  tx„    =    (aw1cosa)1t)x1  (5-10) 

x5  +  TXr    =    (au>«  cosct»jt)Xg  (A-6) 

*(2N+D  +  T*(2N+1)  =  (awl  C0S  ^  X(2N+1)  (A"6a) 
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Each  successive  term  in  the  series  of  equation  (5-12)  can  be  shown  to  be 
smaller  than  the  preceding  term  by  a  factor  of  approximately  a.    (Refer  to 
ref.   10  for  further  discussion.)    Since  a  is  of  the  order  of  2  x  10"    radians, 
this  series  is  seen  to  converge  rapidly.    Hence  all  terms  beyond  x«  can  be 
reasonably  neglected. 

The  following  evaluation  of  the  series  and  associated  coefficients  clearly 
shows  the  appearance  of  the  factor  a  which  permits  terms  beyond  a  certain 
point  to  be  neglected,    (See  Appendix  A  for  solution  details.) 

Aq^qa    =  x  k  =  elastic  restraint  =  0 

x  =  xi  +  x3  +  x5  +  other  terms  which  are  neglected  (5-12) 

x-    =  A  sin  aj«t  +  B  cos  a>2t  (A -7) 

x«  =  C  sin  a>3t  +  D  sin  w.t  +  E  cos  ct?«t  +  F  cos  w^t  (A- 17) 

x5  =  G  sin(2w1t+  a>2t)  +  ^cos^a^t-h  a>2t)  +  J  sin  (2^  -  a>2) 

+  Kcos(2w1t-  u>2t)  +  Lsin(-u>2t)  +  Mcos(-a>2t) 

+  Nsin(w  t)  +  Pcos(a>9t)  (A-37) 

2  z 


b  +  cto>2  c-brct>2 


(A-12)  B  =  -^  (A-14) 


2     2  2    2 

1+TCJ2  l+TOJ2 

2lwAB-tAwJ  2lu)AB  +  tAojJ 

C  =  i ?-        (A -24)  D  =    * * (A-28) 

2a>3(l  +  r2a>2)  2w4(l  +  r2u>2) 

-aw^A  +  rBw,)  ao^A  -  tBw4) 

E  =    t 1_      (A-26)  F  =    * —         (A-30) 

2o3(l  +  T2^)  2w4(l+  r2o>2) 
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G  = 


ao>1  [E  -  rC(2ct)1  +  u^)  ] 
2(2^  +  a>2)  [1  +  t2(2oj1  +  o>2)2  ] 


-aojj  [C  +  rE(2w1  +  w2)  ] 
2(20^  +  w2)  [1  +  r2(2w1  +  oj2)2 


(A-47) 


(A- 48) 


ao>1[F-  rD(2cd1  -o?2)  ] 
2(2o»1  -  a>2)  [1  +  t2(2Wi  -  w2)2  ] 


K  = 


-aw1[D+  tF(2o)1-o)2)] 
2(2^  -  cj2)  [1  +  t2(2w1  -  oj2)2  ] 


(A- 50) 


(A-49) 


L  = 


aw*  [-E  +  rCcu2 
2o>2[l  +  r2o>2] 


(A- 51) 


M 


-aw-  [C  +  TEct>„  ] 
2w2  [1  +  T    a>2  ] 


(A- 52) 


N 


acjj  [F  +  tDw2] 
2a)2[l  +  t2  w2  ] 


(A- 53) 


aw^D  -  tFoj2] 
2w9[l  +  r20Jp] 


(A- 54) 


.0009  sec  =  .  0  x  10     sec 
,-3 


aw,oV   =    2  x  10"     radians 
max 


(f>^ 


'+c2)max=    4  x  10  3  radians 


Typical  values  as  used  by 
the  MIT  Instrumentation 
Laboratory  Angular  Shaker 


With  negligible  elastic  restraint  between  the  float  and  case,  the  series 
solution  consists  of  sinusoids.    Such  sinusoids  will  exhibit  a  time-average 
value  of  zero  except  for  special  cases  which  are  of  interest. 

The  first  case  has  been  studied  in  ref.   10  and  concerns  equal  frequencies, 

wl  =  w2 
As  in  further  explained  in  Appendix  A,  equal  input  frequencies  will  cause  the 
terms  of  D  sin  ojA  and  F  cos  uA  to  have  indeterminate  values.    By  applying 
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L' Hospital's  Rule,  a  constant  drift  term  which  is  time  dependent  is  found 
to  exist.    This  term  is: 


,t  .[j*,] 


Drift    =  |  where  u>  =  w.  =  u>2 


This  existing  drift  is  identical  to  that  found  in  ref.   10. 

The  second  case  of  interest  concerns  the  almost  natural  region  for  sus- 
picion, namely  harmonics.    Terms  concerning  Xp-  of  the  series  give  indetermi 
nate  values  when   2u)^  =  o>2»    By  L'  Hospital's  Rule,  a  drift  term  appears  which 
is  very  similar  to  the  term  arising  from  equal  frequency  input.    This  term  is: 


■m 


Drift   =  1  t  See  Appendix  A  for  proof  and 

necessary  mathematics 

2w«  =  o>2 

If  further  series  terms  such  as  x„  and  Xg  were  evaluated,  drift  terms 
associated  with  other  harmonic  multiples  are  found.    Thus,  the  conclusion  is 
that  equal  frequencies  and  harmonics  are  capable  of  causing  drift  for  the  case 
of  negligible  elastic  restraint. 

The  physical  case  will  probably  display  many  superimposed  input  fre- 
quencies.   This  is  likely  to  have  some  harmonics  thereby  creating  drift.    Such 
drift  is  substantiated  by  this  theoretical  solution  but  it  is  suspected  that  fre- 
quencies approaching  the  harmonics  may  well  create  actual  physical  drift. 

5.  4.  2  Finite  Elastic  Restraint  Between  the  Float  and  Case 
About  the  Output  Axis 

For  the  case  when  k  4  0  ,  it  is  desired  to  find  the  solution  of  the 
following  equation: 

(k-aojj  cosajjt)  x  +  x  +   tx   =   ba^cosaJgt  -  co»2  sinu^t  +  kbsinojgt  H-kccosco^t 

(5-4) 
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The  mathematical  solution  with  details  is  in  Appendix  A.    The  resulting  series 
with  values  of  coefficients  in  terms  of  the  known  constants  is  as  follows: 

A(I-F)OA  =  x  k  /  ° 

x  =  Xj  +  x3  +  x5  +    other  terms  which  are  neglected  (5-12) 

x^  =      A  sin  ojgt  +  B  cos  ct>2t  (A-67) 

x„  =    D  sin  (joA  sin  wJl  +  E  sin  wA  cos  o>2t 

+  F  cos  cjjt  sin  o>2t  +  G  cos  uA  cos  a>2t  (A-75) 

Xg  =     Hj  sin  2u)A  sin  w2t  +  I  sin  2uA  cos  a>2t 
+  J  cos  2cjjt  sin  qjA  +  L  cos  2aKt  cos  u)A 
+  M  sin  ojgt  +  N  cos  cu2t  (A-89) 

Expressions  for  the  coefficients  in  terms  of  known  constants  may  be  found  in 
Appendix  A.     The  coefficients  have  been  left  in  determinant  form. 


By  the  use  of  suitable  trigonometry  relations  the  following  forms  are 
x«    =       —  [cos(wj  -  w2)t]    -   — [cos(w«  +  w2)t] 

Li  Li 


obtained  for  x«: 


E  E 

+  —  [sin(a>.  +a>2)t]     +  [sin(a>..  -  a>2)t] 

Li  Li 

F  F 

+  —   [sin(a^  +  a>2)t]    -   —  [sin(a>j  -  co2)t] 

c  c 

+  —  [cos(a>*+ a>2)t]     +    — [cos(cj- -  cojt]         (A-75a) 

Li  Li 

When  w.,=  o>2,  the  case  of  ref.   10  occurs.    An  indeterminate  form  does  not 
appear  because  k  is  finite.    Thus  time  dependent  drifting  is  not  a  factor  with 
non-zero  elastic  restraint.    However,  this  special  case  of  equal  frequencies 
gives  a  value  of  unity  to  terms  containing  cos(a^  -  a>2)t  and  the  result  is  a 
steady  non-zero  value  for  the  angle  between  the  float  and  inertial  space  about 
the  output  axis.     This  is  verified  by  ref.   10. 
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The  other  case  of  interest  with  finite  elastic  restraint  again  concerns 
harmonics.    This  can  easily  be  seen  by  converting  x5  terms  by  means  of 
trigonometry  relations  to  the  following: 

x,    =        -^L[cos(2Wl-o)9)tI     .   ^1   [cos(2«1  +  w2)t] 
o  2  i        £  2 

+   —    [sin(2^1  +  a>2)tl      +     —    [sin(2oJ1  -  o^Jt] 
+    —   [sin(2w1+w2)t]     -    —    [sin(2o;1  -  w2)t] 

A  2 

+    —   [cos(2w1  +  o)2)t]     +     —    [cos(2w1  -  w2)t] 
2  2 

+     M  sin  u>2t  +   N  cos  o>2t  (A-89a) 

The  drift  resulting  from  the  harmonic  when  2a).  =  cj2  is  clearly  seen  as 
a  possibility.    However,  because  of  a  finite  value  for  elastic  restraint,  indeter- 
minate forms  do  not  appear.    Nevertheless,  this  particular  condition  causes  an 
angle  to  take  place  about  the  output  axis  between  the  float  and  inertial  space. 
This  occurs  when  cos(2o>1  -  ct>2K  is  unity  in  value. 

This  case  gives  a  term  showing  that  a  constant  angular  displacement 
occurs  for  the  harmonic  condition  instead  of  a  time  dependent  drift  term. 

Such  a  condition  is  almost  a  natural  situation  for  a  gyroscope.    As  the 
finite  restraint  is  increased,  the  displacement  between  the  float  and  the  refer- 
ence in  inertial  space  should  become  smaller.     The  extreme  limits  are  for 
zero  elastic  restraint  and  infinite  elastic  restraint.    The  former  case  gives  a 
time  dependent  drift  while  k  =  oo  creates  a  situation  permitting  no  angle,  what- 
soever, to  exist  between  the  float  and  inertial  reference. 

The  most  interest  undoubtedly  lies  with  the  case  of  finite,  non-zero, 
elastic  restraint.    This  is  the  existing  case  of  practical  importance.    The 
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concluding  results  are  that  1)  a  time -independent  drift  term  is  possible  if 
input  frequencies  contain  harmonics  or  equal  values,    2)  float  angular  dis- 
placement with  reference  to  inertial  space  depends  on  the  degree  of  elastic 
restraint  in  the  gyroscope,  and  3)  high  harmonic  terms  may  possibly  be 
neglected  because  of  their  small  value  which  becomes  insignificant  when 
compared  to  manufacturing  tolerances. 

It  is  pointed  out  that  a  drift  may  appear  to  be  present  at  any  one  instant 
of  time  if  a  time-averaged  measurement  is  not  obtained  even  if  harmonics  are 
not  present.     This  could  be  of  noticeable  magnitude  dependent  upon  the  input 
frequencies  and  their  own  absolute  magnitudes. 

It  is  strongly  recommended  that  this  theory  be  applied  to  the  MIT  Instru- 
mentation Laboratory  angular  shaker  in  an  attempt  to  verify  the  results.    Har- 
monic inputs  appear  feasibile  for  the  shaker  by  minor  adaptation  after  sufficient 
data  has  been  collected  concerning  the  case  of  two  inputs  of  equal  frequency. 

Terms  in  the  approximate  series  solution  beyond  Xr  have  only  theoretical 
interest  at  present.     Until  gyroscopic  tolerances  exceed  their  present  values, 
common  sense  tells  us  that  the  selected  series  is  as  good  for  analysis  as  would 
be  the  unknown  closed  form  solution  for  the  equation  of  motion. 

5.  5     Drift  of  a  Typical  Gyroscope  for  the  Case  of  Negligible  Elastic 

Restraint  About  the  Output  Axis  with  Two  Equal  Input  Frequencies 

The  ensuing  results  are  based  on  the  series  solution  for  the  equation  of 
motion  when  the  elastic  restraint  between  the  case  and  float  about  the  output 
axis  is  negligible.    A  further  stipulation  is  that  o>1  =  a^    in  order  that  this  in- 
formation may  apply  directly  to  the  MIT  Instrumentation  Laboratory  angular 
shaker. 
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The  angular  shaker  has  a  frequency  range  from  twenty  (20)  cycles  per 
second  to  one  hundred  (100)  cycles  per  second.  The  constant  "a"  has  a  max- 
imum value  of  .  002  radians  and  the  constants  "b"  and  "c"  are  limited  such 

J~~ 2        JT* 
b    +  c     has  a  maximum  of  .  004  radians.    A  typical  value  of  r,  the 


time  constant  of  the  damper-mass  system,  is  .  0009  seconds. 
By  definition: 


L(I-C)NS 


a  sin  ojA 


A(I-C)OA   =   b  sin  w2t  +  c  cos  U)2t 


(5-2) 
(5-3) 


For  this  particular  cases    cj  =  on  =  coo  as  the  angular  vibrator  is  capable  of 
applying  only  one  frequency  at  a  time  to  the  two  input  axes. 


If  c  =  0,   A/,  r\NS  and  A/,  q\qa   are  in  phase  and  they  will  differ  only 
by  the  magnitudes  of  "a"  and  "b".    (a  =  b  is  a  special  case.)   In  general,  "b" 
and  "e"  will  have  finite  values  thus  causing  A/j  r\j>jg  and  A^  q)qj^  to  differ 
by  a  phase  angle  as  well  as  magnitude.    (Here,  again,    a  =  -J  b2  +  c^    is  a 
special  case.) 


b  sin  art  +  c  cos  art  =|sin(art  +  cot"   b/c)  "Mb    + 
cot"   b/c    =  0 ,  a  phase  angle 
If  o)j  =  o>2,  the  following  is  drift  rate  for  k     0  : 


Drift  Rate  = 


aco) 


1 T  O) 

c 

2     2 

1  +  T      W 


A  plot  is  made  for  drift  rate/ac  for  various  ratios  of  b/c. 


Drift  Rate 
ac 


=     wf 


1-JL     TO,"1 


-         2    2 

1  +   T      OS 


1 TO) 

C 

7~      2    2 

1+    T    W 


(5-14) 


(5-15) 


(5-16) 
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Positive  and  negative  values  of  b/c  were  selected  to  correspond  to  phase  angles 
from  approximately  ±  5°  to  ±  84°„  Ratios  for  b/c  selected  are  indicated  on  the 
graph. 

The  curves  of  Fig.   5-2  show  that  for  a  frequency  range  of  20  to  100  cycles 
per  second,  the  drift  rate  will  undergo  a  change  in  direction  for  b/c  =  +  2.  0 
around  90  cycles  per  second.    The  frequency  required  for  this  reversal  is  fur- 
ther decreased  for  increasing  values  of  b/c. 

Although  the  curves  represent  specific  values  for  constants,  it  is  believed 
to  be  typical  of  gyroscopes  which  may  be  analyzed  by  the  MIT  angular  shaker. 
It  is  suggested  that  these  curves  be  compared  to  actual  data  from  the  vibrator 
test  machine  at  the  earliest  opportunity. 
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CHAPTER  6 


CONCLUSIONS  AND  RECOMMENDATIONS 


6.  1     Conclusions 

1.  When  a  floated  single -degree -of -freedom  integrating  gyro- 
scope is  subjected  to  vibrations  undesired  torques  or  drift 
rates  may  exist.    A  summary  of  such  unwanted  errors, 
arising  when  the  vibrations  are  angular  in  nature,  is  pre- 
sented in  table  form.    The  errors  are  classified  by  name, 
type  and  cause. 


Name 


Anisoinertia  effect 


Tyfre  Error 
Dynamic 


Cross  coupling  coning 


Gyroscopic 


Cause 

Simultaneous  angular 
motions  about  the  case 
input  and  spin  axes 
along  with  unequal  mo- 
ments of  inertia  about 
these  axes. 

Float  precession  about 
its  output  axis  which 
causes  a  component  of 
the  case  spin  axis  angular 
velocity  to  be  coupled  to 
the  float  input  axis  angu- 
lar velocity. 
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Platform  coning  or  Geometric  Simultaneous  angular 

gimbal  walk  velocities  about  plat- 

form axes  such  that 
the  float  input  axis 
undergoes  a  periodic 
conical  motion 

Sculling  effect  Geometric  Simultaneous  angular 

motion  about  any  plat- 
form axis  and  linear 
acceleration  along  a 
perpendicular  axis 
producing  a  rectified 
acceleration  along  an 
axis  perpendicular  to 
both  axes. 

2.  A  second  table  is  presented  to  illustrate  which  effects  are 
present  when  a  single -degree -of -freedom  integrating  gyro- 
scope is  subjected  to  constant  angular  velocities  and  periodic 
vibrations  which  may  be  in  or  out  of  phase. 

Constant  Angular     Periodic  In  Phase      Periodic  Out  of  Phase 
Velocities  Angular  Velocities        Angular  Velocities 

Yes  No 

No  Yes 


No  Yes 


Anisoinertia 

Yes 

Cross  coupling 

Yes 

coning 

Platform  coning 

Yes 

or  gimbal  walk 

Sculling* 

Yes 

Yes  No 


♦For  sculling  the  phase  refers  to  the  phase  relation  between 
the  linear  acceleration  and  the  Euler  angle  motion  (space  angle  motion). 
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3.  Coning  is  a  cause  of  drift  for  a  single-degree -of -freedom 
integrating  gyroscope. 

4.  Cross  coupling  coning  is  inversely  proportional  to  frequency 
and  directly  proportional  to  the  product  of  the  amplitudes  of 
the  angular  motions. 

5.  Cross  coupling  coning  may  be  minimized  by  use  of  a  tight 
servo  loop  which  keeps  the  gimbal  angle  close  to  zero. 

6.  Platform  coning  is  directly  proportional  to  frequency  and  di- 
rectly proportional  to  the  square  of  the  amplitude  of  angular 
motion. 

7.  Sculling  generated  drift  is  directly  proportional  to  the  ampli- 
tude of  the  angular  motion  and  linear  acceleration. 

8.  For  zero  elastic  restraint  between  the  float  and  case  about 
the  output  axis:  (Based  on  a  series  solution  of  a  dynamic  model) 

(a)  Constant  drift  rate  appears  when  two  equal  frequencies 
are  simultaneously  applied  to  the  spin  and  output  axes 
(e.  g. ,    Wj  =  a)2  ) 

(b)  Constant  drift  rate  appears  when  the  frequencies  applied 
simultaneously  to  the  spin  and  output  axes  have  a  har- 
monic relation   (e.g.  ,    2oj^  =  Wg) 

(c)  The  total  drift  at  any  instant  of  time  consists  of  several 
superimposed  sinusoids.    Except  for  the  special  cases 
concerning  inputs  of  equal  or  harmonic  frequencies  re- 
ferred to  in  (a)  and  (b),  the  time  averaged  value  of 
angular  displacement  of  the  float  with  respect  to  inertial 
space  about  the  output  axis  is  zero. 
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9.  For  finite  elastic  restraint  between  the  float  and  case 
about  the  output  axis:  (Based  on  a  series  solution  of  a 
dynamic  model) 

(a)  Drift  terms  which  are  time  independent  appear 
when  two  equal  frequencies  are  simultaneously 
applied  to  the  spin  and  output  axes  (e.  g. ,  a)«  =  ct>2). 
These  drift  terms  have  a  constant  magnitude  de- 
pendent upon  a,  b,  c,  r,  a),,  and  Wn.  but  not 
dependent  upon  time. 

(b)  Drift  terms  which  are  time  independent  appear 
when  the  frequencies  applied  simultaneously  to 
the  spin  and  output  axes  have  a  harmonic  relation 
(e.  go ,    2(0.  =  o>2)°    These  drift  terms  have  a  con- 
stant magnitude  dependent  upon  a,  b,  c,  r,  a**, 
and  a)«  ^u^  n°t  dependent  upon  time. 

(c)  The  total  drift  at  any  instant  of  time  consists  of 
several  superimposed  sinusoids.    Except  for  the 
special  cases  concerning  inputs  of  equal  or  har- 
monic frequencies  referred  to  in  (a)  and  (b),  the 
time  averaged  value  of  angular  displacement  of  the 
float  with  respect  to  inertial  space  about  the  output 
axis  is  zero. 

6. 2      Recommendations 

1.  Verify  theoretical  results  by  use  of  the  MIT  precision  angular 
vibrator.    This  should  be  done  for  various  values  of  elastic 
restraint. 
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2.  The  present  vibrator  is  capable  of  applying  two  simultaneous 
inputs  of  equal  frequency  to  two  orthogonal  gyroscope  axes. 
It  is  suggested  that  the  machine  be  investigated  to  determine 
if  input  frequencies  with  harmonic  relationships  can  be  simul- 
taneously applied. 

3.  Linearity  should  permit  results  obtained  from  a  specific  fre  - 
quency  applied  to  one  axis  to  be  superimposed  upon  the  results 
obtained  from  a  different  frequency  applied  to  another  gyroscope 
axis.    With  th^s  procedure,  it  may  be  possible  to  obtain  physical 
data  concerning  simultaneous  inputs  of  any  desired  frequencies 
within  the  limit  of  the  vibrator.    Information  for  assorted  com  - 
binations  of  input  frequencies  to  two  orthogonal  gyroscope  axes 
could  thus  be  obtained  from  an  angular  vibrator  capable  of 
applying  only  one  frequency  at  a  time. 

4.  Great  care  should  be  taken  to  obtain  valid  data  in  the  low  fre- 
quency range.    This  is  the  suspected  source  of  many  uncertainties 
and  may  be  caused  by  the  relatively  obscure  bending  modes  of 
large  missiles  and  high  speed  aircraft. 
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APPENDIX  A 

SERIES  SOLUTION  FOR  THE  EQUATION  OF  MOTION  OF  A 
SINGLE-DEGREE-OF-FREEDOM  INTEGRATING  GYROSCOPE 


A.  1     Series  Solution  for  Negligible  Elastic  Restraint  Between 
the  Float  and  Case  About  the  Output  Axis   (k-*0) 

The  general  equation  is: 
(k-aa>«  cos  wA)  x  +  x  +  rx  =   bcj«  cos  a^t  -  c a>2  sin  u)Jl  +  kb  sin  cjgt 

+  kc  cos  a>2t    (A-l) 
If  k-*0,  we  have: 

(-aa^  cos  Q)A)  x  +  x  +  t*x  =   ba>2  cos  a>2t  -  ca>2  sin  a>2t  (A-2) 

In  accordance  with  the  series  established  in  Section  5.  5,  the  interest  lies  with 
the  following  terms: 

x  =  xx  +  x3  +  x5  +  .   .  .  .    +  x2m_1  +  x2m  (A-3) 


where: 


Let 


xl  +   TX1  =  ba,2  cos  U)2t~CQt)2  sin  Ct,2t  (A"4) 

x3  +  tx3  =  (aa^  cos  Wjt)Xj  (A-5) 

x5  +   7X5  =  ^awl  C0S  wl^  x3  (A"6) 

X(2N+1)  +  TX(2N+1)    =  (aa,l  cos  wlt)x(2N+l)  (A"6a) 


Xj    =   A  sin  ojgt  +  B  cos  Wot  (A-7) 

A  a>0  cos  0>ot  -  B  w0  sin  u>0t  (A- 7a) 


x 


-         =      £\    Wn    lUi3    We)**  —  *->    Wn    OJLli    Wn 


X 


j     =    ^A  o>|  sin  Wot  -  B  oji  cos  w«t  (A- 7b) 


Substitute  (A-7),  (A-7a),  and  (A-7b)  into  (A-4) 

2  2 

[Aw,  cos  ojgt-  Bcjg  sin  Wgt]   +    r  [-Aw,  sin  w2t  -  Bwg  cos  oj^t' 

=   ba>2  cos  w«t    -  ca^  sin  Wgt 
Equate  coefficients  of  like  terms  and  solve  for  A  and  B. 


2    2 
A  +  7    a>2  A   =   b  +  crojg 


(A-8) 


9 

Aw„  -rBwn    =   bo>2            (from  cosine  terms) 

(A-9) 

2 
Ba>2  -tAu>2    =    -  cu> g        (from  sine  terms) 

(A-10) 

A  =  b  +  r  B  o>2 

(A- 9a) 

B  =  c  -  t  A  a>2 

(A- 10a) 

A  =  b  +  tw„  [c  -  tAoj9] 

(A- 11) 

b  +  cto>2 
"         2~2~ 

1  +    T      W« 


B  =  c-rco2[b+  rBojg] 

2    2 

B  +  t    w0B    =    c  -brcj. 


(A- 12) 


(A- 13) 


B  = 


c  -  br  a>2 
"         2~~2 

1  +   T      0^2 


Only  those  terms  affecting  the  steady  state  solution  are  considered. 
An  equation  for  x«  is  obtained  by  substituting  (A-7)  into  (A- 5). 
x„  +  tx„    =   (aw j  cos  o>1t)(A  sin  Wgt  +  B  cos  o^t) 

x«  +  tx3    =   A  aw.  cos  co A  sin  w2t  +  Bawj  cos  Wjt  cos  o>2t 


(A- 14) 


(A- 15) 
(A- 15a) 
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From  basic  trigonometry  relations,  the  following  are  obtained: 
cos  ajjt  sin  o>2t     =  —   [sin^.t  +  u^t)  -  sin(u)..t  -  w2t)  ] 
cos  Wjt  cos  cjgt     =  —   [cosCcjjt  +  (j) J)  +  cos(a>.,  t  -  Ct)2t)  ] 

Let   w3t  =  (ajjt  +  o>2t)      and   wA  =  {uA  -  a>2t)  t  . 

In  view  of  these  relations: 


Aaaj.,  Baco 


t  M  Aaojj  Baco.. 

x3  +  tx3    = [sin  ojgt  -  sin  o^t]   +  [cos  uA  +  cos  vA]    (A- 16) 


Let 


x„  =   C  sin  u)A  +  D  sin  cu^t  +  E  cos  uA  +  F  cos  uA  (A- 17) 

Xg  =   Cco^  cos  o>gt  +  Dw.  cos  wA  -  Ea)«  sin  Wqt-  Fw.  sin  ojA  (A-17a) 

••2  2  2  2 

x«  =  "^3  sin  (jjA-Dco^.  sinaj.t-Ew^  cos  o^t-Fcj,  cos  w.t  (A- 17b) 

Substitute  (A-17),  (A-17a),  and  (A-17b)  into  (A-16)  . 

[Ccj«  cos  u)A  +  Da>4  cos  a>,t  -  Eci>«  sin  o)«t  -  Fw,  sin  u)A] 

2  2  2  2 

+  [-tCc^3  sin  wA  -tDw,  sin  cj.t  -tEc^  cos  cjgt  -tFw,  cos  <j)A] 

Aslu).  Aslu).  Baco.  Baa). 

=  sin  qjA  - sin  ooA  +  cos  uA  + cos  u)A      (A-18) 

2  6  2  *  2  6  2  q 


Equate  coefficients  of  like  terms  and  solve  for  C,  D,  E,  and  F. 

o 
)A  :         Cct>«  -rEojg 

cos  w4t :  Dcj4  -tFcj^    =     —  (A-20) 


2           Bao;1 
coso)«t:         Cw«  -tEo>«    =     —  (A-19) 


Bawl 


2 


Aaw. 


p  /lacy* 

sinw3t:        -Ew3  -tCw3    =     —  (A-21) 
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Aaw, 


sin  uA  :       -Fw,  -rDo/ 


(A-22) 


C    = 


1/2  Bawj  +  rEo>3 


(A- 19a) 


E    = 


1/2  Aaojj  -  7CW3 


(A-21a) 


Substitute  (A-21a)  into  (A-19) 


Coj0  -tw! 


-1/2  Aawj  -tCo>3 


U)t 


Baw. 


(A-23) 


2-, .3 


Ccj3  ■+   l/2rAao)1w3  +  t^Co>3    =    1/2  Bawj 


2    3 


C(oj3  +  r    w3)    =    1/2  Bao^  -  1/2  rAaw.ojg 


Bacu. 


C    = 


rAaojj  a>3 


2(cu3+  r2^) 


(A-24) 


Substitute  (A- 19a)  into  (A-21) 


-Eoj3  -   ru>3 


1/2  Baojj  +   tEo>3 


Aaw, 


2- ,.3 


-Eo>3  -  1/2  rBaojj  u>3  -tzEo>3    =    1/2  Aac^ 


2    .3 


E(o)3  +  r    w3)    =    -  1/2  Aacjj    -  1/2  rBao^a^ 


(A-25) 


E    = 


-  Aaw,    -    rBaoj.Wg 

2     3 
2(w3  +   r    o>3) 


(A-26) 
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1/2  Baojj  +  tFu>4 


u> 


(A-20a) 


1/2  Aawj  -  rDc4)^ 


w. 


Substitute  (A-22a)  into  (A-20) 


Do>4  -   t  a>4 


2 

1/2  Aacj.  -  tDu)m 


Bacu. 


Do>4  -  1/2  TAao)1o>4  +  r  DcoJ    =    1/2 Baa^ 
D(a>4  +  r    w4)    =    1/2  Baw.   +   1/2  tA&u.u)* 


D    = 


Baw*  +  rAaa^a>4 

2    3 
2(w4  +  r    a>4) 


Substitute  (A- 20a)  into  (A -22) 


■  Fo>4  -   r  w4 


l/2Baw1  +  7Fa>4 


Wi 


F<x>4  -  1/2  rBacu1a>4  -  r   F<x>4 
2.3 


=   -  1/2  Aao) 


1/2  Aaojj 


F(w4  +  t^w4)    =    1/2  Aao)1  -  1/2  rBaa^u^ 


Aacu1  -  TBao)jO)4 
2(a>4  +  r2^) 


(A-22a) 


(A-27) 


(A-28) 


(A  29) 


(A- 30) 


An  equation  for  x5  is  obtained  by  substituting  (A- 17)  into  the  following: 

x5  +  Tx5   =   (aw.  cos  Wj)x3  (A-31) 

Xc  +  TXr    =  aoJ1  cos  w..t[C  sin  wJL  +  D  sin  a>4t  +  E  cos  w^t  +  F  cos  a^t]     (A-32) 


(A- 33a) 

(A-33b) 

(A-33c) 

(A- 33d) 

(A-34) 

(A-35) 

TX5  +  x5   =  Cacui  cos  wi*  sin  Wgt  +  Da.u*  cos  coA  sin  w-t  +  EauK  cos  w.t  cos  a>gt 

+  Fao^  cos  cjjt  cos  coA  (A-32a) 

Trigonometry  Relations 
cos  coA  sin  coA    =  —   [sinC^jt  +  coA)  -  sin(u>A  -  coA)  ] 

cos  (*) A  sin  coA   =  —  [sinCojjt  +  coA)  -  sin{coA  -  coA)  ] 

cos  co A  cos  coA  =  —   fcos(wj, t  +  coA)  +  cos(ok  t  -  coA)  ] 

cos  it) A  cos  co  A  =  —  [cosCaJjt  +  w.t)  +  cos(a>..t  -  coA)  ] 

co  A   -   (coA  +  coA) 

co  A   =    (coA-co2t) 

—        •  i 

TX,  +  x5    =  —  Ca.co*  [sin(2w1t  +  coA)  -  sin(-o>2t)  ] 

+  —  Daaj.  [sin(2w<t  -coA)  -  sin(o>2t)  ] 

+  —  Eacu,  [cos(2coA  +  u)A)  +  cos(-o)2t)  ] 
z 

+  —   Fao^  [cos^ajjt  -o)2t)  +  cos(o)2t)  ]  (A-36) 

Let 
x5    =   G  sin(2a>jt  +  coA)  +  H*  cos^ajjt  +  o>2t)  +  J  sin(2coA  -  co^)  +  K  cos(2aJ1t  -  a>2t) 

+  L  sin(-a>2t)  +  M  cos(-a>2t)  +  N  sin(a>2t)  +  P  cos(a>2t)  (A-37) 

x5   =   0(20;-  +  a>2)  cos(2w1t  +  w2t)  -  H1(2a>1  +  a>2)  sin(2a>jt  +  coA) 
+   J(2a)1  -  cj2)  cos(2a>1t  -  a>2t)  -  K(2a)..  -  uO  sin(2w-t  -  u>2t) 
-    La;2  cos(-a>2t)  +  Ma)„  sin(-w2t)  +  Nw2  cos(a>2t)  -  Pa,'2  sin(a>2t)         (A-37a) 
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••  2  9 

x5    =  -G(2o>1  +  a>2)    sin(2wjt  +  w2t)  -H^a^  +  a>2)    cos^a^t  +  a>2t) 

-  J(2w1  -  cd2)    sin^a^t  -  a>2t)  -  K^Wj  -  a>2)2  cos^a^t  -  c^2t) 

2  2  9  9 

-Lo)2  sin(-a>2t)  -  M  a>2  cos(-a>2t)  -Na>2  sin(a>2t)  -Pa>2  cos(aJ2t) 

(A-37b) 

Substitute  (A-37),  (A-37a)  and  (A-37b)  into  (A-36). 

2  2 

-rG(2o)1  +  w2)    sin(2w1t  +  a>2t)  -t  H1(2co1  +  a>2)    cos^a^t  +.  ^2t) 

2  2 

^rJ(2a)1  -  a>2)    sin(2a>jt-  a>2t)  -  t K(2a>..  -  a>2)    cos(2co<t  -  u>2t) 

2  2  2  2 

-rLa)2sin(-a)2t)   -  tM  aJ2cos(-Ci>2t)  -r  Nw«  sin(a>2t)  -  r  P  a>2  cos(aJ2t) 

+    G(2ok  +  oj2)  cos(2ct>«t  +  a)2t)  -  H^a*.  +  a>2)  sin(2ct>1t  +  aj2t) 

+    J(2cjj  -  a>2)  cos(2ct)^t-  a>2t)  -  K(2w.  -  a>2)  sin  (2a). t-  a>2t)  -La>2  cos(-a>2t) 

+  Ma>2  sin(-a)2t)  +  Na>2  cos(co2t)  -  Pa)2  sin(a>2t)     =     1/2  Caa^  sin(2o>1  +  a>2t) 

-  1/2  Cawj  sin(-w2t)  +  1/2  Daw.  sin(2a>..  t -out)  -  1/2  Daa^  sin(a>2t) 

+  1/2  Eaa)-  cos(2a>1t+  a>2t)  +  1/2  Eaaj.  cos(-cj2t)  +  1/2  Fact),  cos^oj.t -a>2t) 

+  1/2  Faco1  cos(a>2t)  (A- 38) 

Now  equate  like  parts  and  solve  for  coefficients. 

sin(2w-t+  o>2t):       -rG(2w1  +  cu2)2  -  H1(2a>1  + o>2)    =  —  Cau^  (A-39) 

z 

cos(2w1t+  w2t):       -tH1(2o)1  +  cu2)2  +  G(2o;1  +  o>2)   =  —  Eaa^  (A-40) 

sin(2w1t-w2t):        -7,1(20^  -  a>2)2  -K^ojj  -  a>2)     =   J_  Daa^  (A-41) 

z 

cos^o^t-o^t):       -tK^o^-u^)2  +  JC2WJ  -  a>2)    =   —  Fac^  (A-42) 

sin(-0J2t):  -tLw2  +  Mo>2    =    -  1/2  Caw x  (A-43) 


68 


cos(-u)2t) 

sin  (o>2t)  : 
cos  (cjgt)  • 


tMu,  -Loj«    =  — Eaa), 


2  1 

•TNcdg-P^o    =  —  Daa>< 
2 

2  1 

•  tPo>2  +  Na>2    =    —  Faa^ 
2 


-  1/2  Caw*  -  r  G^u*  +  oj0)' 
Ht    -    i 1 — 

Substitute  (A- 39a)  into  (A -40)  and  solve  for  G: 

,2 


(A-44) 

(A45) 
(A-46) 

(A- 39a) 


■  t  (2ct>..  +  w2)' 


-  1/2  Caa^  -  r  G(2w1  +  w2)' 
(2o>1  +  w2) 


+  G(2o)1  +  a>2)  = —  Eawj 


(A- 40) 


1  2  3  1 

—  rCacoj  (20^  +  Wg)  +  t-  G(2oj1  +  w^)    +  G(2o)1  +  cj2)  =  —  Eaajj 

2  2 

G     (20^  +  a)2)  +  t2(2u>1  +  u>2)3       =  —  Eact>1  -—  TCaw1(2a>1  +  a>2) 
Eaw.  -  rCacj^  (2oj^  +  u)^ 


2  3 

2[(2w1  +  o>2)  +  T    (2wj  +  ol>2)    ] 


(A-47) 


or 

G    = 

awjE—  rC(2w1  +  ci>2)  ] 

2[(2a>1  +  a>2)  +  t2(2w1  +  a>2)3] 

G    = 


1/2  Eaa^  +  t  H1(2u>1  +  a>2)' 


(A-47) 


(A -40a) 


-t{2(j).  +  o)2)' 


(2d).   +  Wg) 

Into  (A- 39) 

l/2Eao>1  +  t  H1(2w1  +  o>2V 


Substitute  (A~40a)  into  (A- 39)  and  solve  for  Hj  : 

,2 


(20^  +  0^2) 


(A-39) 
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-  —  rEao^  {2u1  +  u>2)  -  t  H1(2a>1  +  u>2)3  -  H^u^  +  a>2)    =  -i  Cawj 
^  2 


H1[(2w1  +  a>2)  +  t2(2o)1  +  cj2)3] 


—  Cao), r  Eaa)1(2a>1  +  a>2) 


-Cawj  -r  Ea.u)A2cj^  +  o>2) 
2[(2o>1  +  cj2)  +  t2(2o;1  +  w2)3 


(A-48) 


or 

Hl 

aojj  [-  C  -r  £(2^).  +  a>2)  ] 

2[(2o>1  +  a>2)  +  r2(2w1  +  w2)3] 

J  = 

l/2Fact>1+  tK(2oj1-w2)2 

(2^  -  o>2) 


(A-48) 


(A-42a) 


Substitute  (A- 42a)  into  (A- 41)  and  solve  for  K : 

2 


7(205.   -  0>2)' 


1/2¥slu)1  +  tK(2u)1-u)2)' 
{2w1  -  o>2) 


-  K(2o>1  -  cj2)  =  —  Dac^    (A-41) 


1  2  3  1 

—   TFaojj  (2o>1 -o)2)  -  r   K(2oJ1-ct))     -  K^o^  -  o>2)    =—  Dao>1 


K[  (2o>1  -  w2)  +  r2(2o)1  -  o>2)3  ]    =  -  i.  Dao^  -  —  T  Fao)1(2w1  -  u>2) 


K    = 


aojj  [-D-r  F(2w1  -o>2)  ] 
2[(2o»1  -  u)2)  +  r2(2u)i  -  o>2)3  ] 


(A -49) 


-  1/2  D3wx  -   tJ(2u1  -  u>2)' 
(2o)1  -  Q)2) 


(A-41a) 
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Substitute  (A-41a)  into  (A-42)  and  solve  for  J 


t(2oj1  -  u>2)' 


-  1/2  Daw1  -  r  J{2u1  -  u)2)' 


+  JC20).  -  oj2)    =  —  Facjj 


(2^  -  u)2) 
—  r  Dao)  (20^  -  a>2)  +   r  J  (2a?1  -  a>2)3  +  J(2c^1  -  a>2)  =  —  Facoj 


2  3  1  1 

J[(2a>1-w2)+   T   (2cu1-oj2)    ]    =  —  Fawj  -  —  rDau^  (20^  -  a>2) 


J        = 


acjj  [F  -  tD(2w1  -  cj2)  J 
2[  (2o>1  -  cu2)  +  r  2(2w1  -  a>2)3  ] 


M 


1/2  Caw.  +  rhu>2 


(A-42) 


(A- 50) 


(A-43a) 


U)c 


Substitute  (A-43a)  into  (A-44)  and  solve  for  L  : 


-T  Ct>< 


v 

1/2  Caw..  +  tLcj2 


Lw0  =  —  Eaa).. 
Z      2  X 


1  2        3  1 

—  TCaoKCJg-r     Lw2   -  Lo>2    = —  Eao^ 

2  ^ 

03  i  i 

L(o>2  +  t *o>2)    =    -  —  Eacj1  +  —  r  Caw^g 

z  z 


L 

awJ-E  +  t  Cct>2] 

2[o,2  +  r2a,|] 

T, 

-l^Eaaij-TMcog 

(A-44) 


(A-51) 


(A- 44a) 
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Substitute  (A-44a)  into  (A-43)  and  solve  for  M  : 


■T  CO, 


1/2  Eao>1  -tMw! 


w« 


Moj< 


Caw. 


-    r  EawjCjg  +   t    Mwq  +  Mw^    =    - —  Cacu* 
2  2 


(A-43) 


M(a>2  +   t   o)2) 


Caw- 


r  Eaa>..a> 


lw2 


M    = 

acu«  (-C  -  r  Ea^) 

2    3 
2(a>2  +  t    a^) 

P    = 

l/^Dawj-  tNo>2 

(A- 52) 


(A-45a) 


Substitute  (A-45a)  into  (A-46)  and  solve  for  N: 


T  0>< 


l^Daa^-TNa^ 


a). 


Ncu, 


Faco. 


(A-46) 


1  2       3 

t  Dao^cjg  +  t  Nwo  +  No>2    =  — 

2  2 


Fact). 


N(o>2  +  t2  a>2 )    =  —  Faa^  +-   tD 
2  2 


lw2 


N    = 


aa>1  (F  +  tDcu2) 

2     3 
2(a>2  +  t    ^2) 


(A- 53) 


N 


1/2  Faajj  +  tPw2 


(A-46a) 
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Substitute  (A-46a)  into  (A-45)  and  solve  for  P  : 


1/2  Fawx  +  r  P  u)\ 


-p"2    = 


Daw, 


(A-45) 


2        3  1 

—  t  Facu^  -  t    Pa>2    -  Po>2    =    -—  Daa^ 
^  2 


P(W9+   T2U)b 


Dact)1 r  FacJiW 


lw2 


awj  (D  -  r  Fa>2) 
2(w2  +  t2co^) 


x  =  xx  +  x3  +  x5 


^m-l  +  *2] 


Xj  =    A  sin  a>2t  +  B  cos  wJ. 


x«  =   C  sin  w«t  +  D  sin  wA  +  E  cos  uJ.  +  F  cos  coA 


(A- 54) 

(A-3) 
(A- 7) 

(A-17) 


x5  =     G  sin(2a)1t+  a>2t)  +  Hj  cos(2u>1t+  a>2t)  +  Jsin(2a);.t  -  a>2t)  +  Kcos  (2a>..t  -  wJ 
+    Lsin(-a)2t)  +  Mcos(-w2t)  +  Nsin(a>2t)  +  Pcos(a>2t)  (A- 37) 

Terms  beyond  Xr  may  be  neglected  if  the  order  of  magnitude  of  a  is  much 
smaller  than  unity  (maximum  a  on  the  MIT  angular  shaker  is  2  x  10" 

radian).    If  x.,  is  defined  to  be  of  the  order  of  magnitude  of  one,  x,  is  of  the 

o 
order  of  magnitude  of  a,  and  Xr  is  of  the  order  of  magnitude  of  a  .    Hence 

each  successive  term  in  the  series  of  x  is  smaller  than  its  predecessor  in 

accordance  with  the  value  of  a. 

For  further  discussion  on  convergence  and  a  process  for  iteration, 
ref.   10  is  recommended. 
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If  only  the  terms  associated  with  x-,  x«,  and  x-  are  considered,  it  is 
readily  apparent  that  the  time  average  of  drift  for  the  sinusoids  is  zero.   Close 
scrutiny,  however,  shows  drift  possibilities  as   u*  and  u>2  approach  one  another 
in  value  and  also  when  harmonics  are  present. 

The  specific  case  where  w.  =  w„  has  been  investigated  in  ref.  10.    The 
general  case  which  is  the  subject  of  this  chapter  degenerates  to  the  specific 
case  when  cu-  =  a>2.    Under  such  circumstances  D  sin  uA  and  F  cos  wA 
require  L'  Hospital's  rule  for  proper  evaluation,  but  eventually  yield  a  drift 
term  which  is  a  direct  function  of  time. 

The  other  case  of  interest  concerns  harmonics.  The  series  terms  asso- 
ciated with  Xe  require  L'  Hospital's  rule  for  evaluation  when  2u)<  =  a>2.  Under 
such  a  condition  drift  terms  appear  which  are  direct  functions  of  time. 

An  example  of  drift  caused  by  the  existence  of  harmonics  such  as  when 
2o)i  =w«  is  shown  for  the  term  J  sin  (2a>«t  -  cjgt)  . 

If   2WJ-*  a>2,  then  (20^  -  a>2)  "*   0  (A-55) 

au>«  [F-tD(2oj1  -  a>9)  ] 
J  sin  (20)^  -  a>2t)    =    i * 1 sin  (2a>jt  -  a>2t)        (A-  56) 

2(2w1-w2)[l+  T2(2Wl-a>2)2] 
Let   (2w1t-w2t)    =   (2w1-a)2)t   =   yt  (A-57) 

(A- 58) 


lim     J  sin  yt   =  — 
y*0                        0 

apply  L'  Hospital's  Rule 

—  [aa^F-  rDy)  sin  yt] 
dy 

ao). 
2 

-  r  D  sin  yt   +  (F  -  rDy)(t  cos  yt) 

-A.  [2y(l  +  r2y2)] 

(l+r2y2)   ♦  y(2r2y) 

(A- 59) 
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aw.. 

lim     J  sin  yt    =   — 

y-*0  2 


(F)(t) 


1 


FacoA 

L_  (A -60) 


If  terms  for  x„,  xg,  etc.  were  evaluated,  drift  terms  would  appear  for  higher 
harmonics.    As  the  series  has  an  infinite  number  of  terms,  it  seems  apparent 
that  any  two  given  input  frequencies  might  create  some  type  of  harmonic  drift , 
but  the  magnitude  of  "a"  associated  with  such  drift  may  make  the  consideration 
of  the  value  a  folly  when  compared  to  the  manufacturing  tolerances  used  when 
constructing  gyros. 

The  case  of  negligible  elastic  restraint  shows  the  possibility  of  drift  terms. 
However,  for  such  a  gyroscope,  the  shortened  series  representing  drift  is  valid 
only  when  values  of  a,  b,  and  c  are  small  compared  to  unity.     Furthermore, 
man  can  only  produce  a  gyro  which  is  not  perfect.    Hence  terms  in  the  series 
which  are  of  less  consequence  than  the  manufacturing  tolerances  will  have  only 
a  theoretical  importance  and  should  not  enter  practical  considerations. 

A.  2     Series  Solution  for  Finite  Elastic  Restraint  Between 
the  Float  and  Case  About  the  Output  Axis   (k  4  0) 

The  general  equation  is  : 

(k  -  ao>1  cos  wA)  x  +  x .+  rx  =   ba^  cos  ^o*  ~  CCt,2  sin  ^2*  +  ^  sin  ^2* 

+  kc  cos  out     (A-l) 

In  accordance  with  the  series  established  in  Section  5.  6,  the  interest  lies  with 
the  following  terms  : 


X    =    Xi     +    X0    +    X 


•  ^m-l  +  x2m  (A-3) 


where: 


1  T  A3  T  ~5  T  *  •  •  "  ~2m 

kx.   +  x-  +  t'x<    =   bo>2  cos  wA  -  ccj2  sin  Wgt  +  kb  sin  wA  +  kc  cos  wA 

(A-61) 
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•  •• 

kx3  +  x„  +  tx,    =   (aa).  cos  ojA)  x.  (A-62) 

kxg  +  Xp-  +  Tx5    =   (aw.,  cos  a)A)x^  (A-63) 

*  •• 

(k-  awj  cos  w^t)  x  +  x  +  t  x   =    bw«  cos  w2t  -  cw«  sin  aj«t  +  kb  sin  u>2t 

+  kc  cos  wA   (A-l) 

Let: 

x  =  x..  +  x«  (A-64) 

where: 

•  •• 

kx.   +  x.   +  TX-    =   bw2cosa>2t  -  ca^sina^t  +  kbsina>2t  +  kccosa>2t 

(A-65) 
Substitute  (A-64)  and  (A-65)  into  (A-l). 

•        •  ••       •• 

(k-aoj.  cosw-t)  (x-  +  Xg)  +  (x..  +  x2)  +  t  (xj  +  x,)    =  bcu^cosa^t  -ca>2sinw2t 

+  kbsinojgt  +  kccoso^t  (A-66) 


or 


•  •  • 

(k  -  aw*  cos  Q)A)  x,  +  x«  +  tx«  =  (aw.  cos  uA)  x^  (A-66a) 

Try  a  solution  for  (A-64)  of : 

x.    =   A  sinojgt  +  B  cos  u)A  (A-67) 

x.    =   AWrtCOS  u)A  -  Ba>2sina>2t  (A- 67a) 

••2  2 

x..    =   Aw,  sincjgt  -Bw2  cosa^t  (A-67b) 

Substitute  (A-67a)  and  (A-67b)  into  (A-65) . 

o 

k  A  sin  w«t  +  k  B  cos  u)Jt  +  Aw,  cos  a>2t  -Bw2  sin  a^t  -  T  A  ««  sin  w2t 

o 
-  tBcJqCOS^^  =   ba>2cosw2t  -  cw2sinw2t  +  kbsinw2t  +  kccosco^t      (A-65) 

Equate  coefficients  of  like  terms. 

2 
sin  ct>2t :  kA  -  Ba>2  -  tAo>2  =    -  ca>2  +  kb  (A- 68) 

2 
cos  a>2t :  kB  +  Aa>2  -  tBo>2  =    bu>2  +  kc  (A- 69) 
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A(k-ro>9)    =   Bo>0  -cct>9  +  kb 


A   = 


rBa>2  -kB  +  ba>2  +  kc 


w« 


rBa>, 


2    "w2 


kB 


it). 


+  b  + 


kc 


(A- 68a) 
(A-69a) 


Substitute  (A- 69a)  into  (A -68 a)  and  multiply  through  by  a>9  : 

kTBa^  -k2B  +  kba>2  +  k2C  -  t2&u\  +  rkBa)2  -  rb(x>\  -  rkcu)2 

2  2 

Boj„  -  coJq  +  kba)« 


(A-70) 


B[    tw2  -k     -r    w«  +  T.ko)«  -w«]    =    rbwp  - k    c  +  r  kccc)p  -ecu! 


B   = 


rbw0  -C/  — A      +  7 

I  "2  J 


kc  -  c 


+  /jL\1_rft,chL+JL 


^ 


i; 


C  it), 


2k 


(:k\2       2    2, 


%)--■ 


c 
B    - 

"i+  jl]2  _T(02crA+jLi 

w2j                   Lc      "8  J 

'*[(;)--] 

2 

(A-71) 


B 


Ak  -  A  t  w«  +  co)„  -  kb 


U), 


Ak 


c*>» 


Arajp+  c 


'2  "2 

Substitute  (A- 68b)  into  (A-69a)  and  solve  for  A  : 


-i^_  (A-68b) 


(TU)Z   -k) 


w. 


Ak        A  „      kb 

Wo  Wo 


+  b  + 


kc 


(A -69) 
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A    =    TAk-  At2ojo+  Tca>0-Tkb-^L_  +  Akr  -  —  +  JLL_  +  b  +  *£- 


1  -rk 


+  r24  +fjL\2.rk 


=  r  cct>2  "T  kb  + 


k'b 


Ct>< 


+  b 


i=f(t)!WHfe)l 

'*[(i)—=]2 

(k-acj.  coso).t)x2  +  x«  +  t  x«    =   (ao^  cosco^t)  x^ 


Let: 
where: 


x2  -  x3  +  x4 


kx3  +  x,+  tx„  =  (aojj  cos  a^t)  Xj 


(A- 72) 

(A- 66a) 
(A-73) 

(A-74) 


Substitute  (A-73)  and  (A-74)  into  (A-66a). 
(k  -  aujj  cos  o>1t)(x3  +  x4)  +  (xg  +  x4)  +  r  (x3  +  x4)    =   (aojj  cos  a^t)  Xj        (A-66a) 


or 


•  •• 


(A-66a) 


(k  -  acuj  cos  (j)  J)  x4+x.+  rx4    =   (aa^  cos  a^t)  x3 

A  solution  is  desired  for  the  following  equation: 

kx„  +  x„+  tx„    =    aw1Acoscu1tsina>2t  +  au>1Bcosa)1tcos  w2t  (A-74a) 

Let: 

Xo  =  Dsincu^tsinoJnt  +  E  sinw.tcos  cjgt 


+  Fcos  cj-itsina^t  +  Gcosa^tcosa^l 


(A-75) 
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Xg    =  Dojj  cosw^tsinWrtt  +  Dw«  sinwjtcosojgt 

+  Eoj^  cos  a),  t  cos  (j)Jl  -EaJn  sin  ooA  sin  ooA 

~  Fa^  sin  Wjt  sin a)2t  +  Fa>2  coscu.tcos  cjgt 

-  Ga^  sincjjtcoscjgt  -  Gw2  cosaKtsinaut  (A-75a) 

••  2 

x3  =  -Dwj  sinco-tsinwgt  +  Tioo+oo^  cosoj.tcoso^t 

2 
+Dw^2  cos  w^t  cos  ooA-  Dc^2  sin  o^t  sin  u)Jl 

o 

-Eo^  sin  ctj^t  cos  w2t  -  Eo^cjg  cosw.tsinojgt 

2 
-E oo .  ojq  cos  Wj t  sin  oo<A  -  Ea>2  sin  a>..  t  cos  a>2t 

2 
-Fa)j  cos  ooAsmooA  -  Fa)^ a>2  sin  cujtcos.Unt 

2 

-Fu>..a>2  sina,itcosa,2t  ~  Fw2  cosa,ltsinw2t 

2 
-Gajj  coscjjtcos  w«t  +  Gcj.  c^2  sinco^tsinw^t 

2 

+  Ga^a^  sinw^tsinajgt  -  Gw„  cosw-tcosajgt  (A-75b) 

Now  substitute  (A-75),  (A-75a)  and  (A-75b)  into  (A-74a)  and  solve  for  D,  E, 
F,  and  G. 

(kD  sin  ooA  sin  wA  +  kE  sin  ooA  cos  ooA  +  kF  cos  ooA.  sin  cjgt  +  kG  cos  ooA.  cos  cjgt) 

+  (Dov.,  cos  oo A  sin  a^t  +  Da^  sin  ooA.  cos  a^t  +  Ect^  cos  cj.  t  cos  u^t  -Ew,  sin  a^ t  sin  a^t 

-Fa**  sin  oj1  t  sin  ooA  +  Fw«  cos  u)A  cos  cugt  -Gw.  sin  ojA.  cos  w«t  -£^2  cos  wl*  s*n  w2^ 

2 

+  ( -rD  Wj  sinw^tsinwgt  +  rDa;..^  cosw-tcosa^t  +  rDa)^2  cos  witcos  w2t 

2  2 

-TDcOgSinaj-jtsin  a^t  -  rE  cj-  sinaj..tcos  c^t  -  rEw.Wn  cos  ajjtsinajgt 

2  2 

-rEcOjWgCOsaj-tsincUgt  -  rE  ajg  sina^tcosa^t  -  tFoj.  cosa>«t  sin  Wgt 

2 
-TFa^a^sinco-tcosu^t  -  rFcj.au  sin  00+  t  cos  a^t  -  tF  a^  cosa^t  sinajgt 

2 

— rGw -  cos  qo A  cos  aj«t+   rGcjjWg  sincj.tsinwgt  +   rGcjjWg  sin  ajjt  sin  out 

2 
-TGcjgCOscJwtcoscjgt)    =   Aac^j  coscj^tsinc^gt  +  Baw.coscu.tcosojgt       (A-76) 

79 


Now  equate  like  terms  of  the  various  trigonometric  function  combinations 
and  solve  for  D,  E,  F,  and  G.    (Note:   A  and  B  are  known) 


(sinojjtsincjgt) 


2  2 

kD  -  Eo>2  -  Fajj  -  tD  w«  -  rD  cj«  +   TGaj-aj^  +   rGa).,  a>2  =  0 

kD  -Ea)2  -Fa>1  -  D(r  cu^  +   r  u>2)  +  2x0^^2  =  0 

D(k-  rwj-rwg)  -Ecjg-Fajj  +  G(2ra)1w2)    =0  (A-77) 


(sincj*tcosa>2t) 


kE  +  Dw2  -  Gcu.  -tEo),-  tE  a>2  -  rFu}..a>2  -  TFa>«a>2 


2  2 

Da>2  +  E(k-  t  co1  -  r  cj2)  -  F^rWjWg)  -  Ga^  =  0  (A-78) 


(cos  cjjt  sin  a>2t) 


2  2 

kF  +  Doj- -Goj2  -  t  Ecu-.ajg  -  rEw-Wg  -  tFcu-  -  tFo)2    =   Aaa^ 

2  2 

Dojj  -  E(2  rcu.a^)  +  F(k  -  r  cjJ  -  t  o>2)  -  Ga>2  =  Aawj  (A-79) 


(cos  u)A  cos  o)2t) 


2  2 

kG  +  Ecu.  +  Fcj2  +  rDco-cUg  +   tDco..u>2  -  tG  o^  -  rG  u>2  =  Baa^ 

D(2to)1cd2)  +  Ecu1  +  Fa>2  +  G(k  -  r  coj  -  t  a>2)   =  Bacuj  (A-80) 


D(k-TO)pTW2)-  Eco2  -  Fcuj  +  G(2r o^o^)  =  0  (A-77) 

Da>2  +  E(k-  to)|  -t  u>2)   -  F(2tcu1w2)  -  Gco1  =  0  (A-78) 

2  2 

Dcjj  +  E^TC^ct^)  +  F(k-r  co^-t  co2)   -  Gco2  =  Aawj  (A-79) 

0 
D(2ro)1ct)2)  +  Eojj  +  Fo)2  +  G(k-r  o>1)  =  Baa^  (A-80) 

Solve  for  D,  E,   F,  and  G  by  determinants. 
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D  = 


0 

0 
Aao) 
Bacj 


2  2 

(k-    TCJ1  -TCt>2) 


270)^2 


Ci>, 


■w. 


2rw^ct)2 


(A-81) 
27  0)<o)2 


■w. 


(k-  T^-  70)2) 


Cl>, 


2  2 

(k-  7  0).  -  7a>2) 


(k-TWJ-T«*) 


CO. 


O). 


O), 


2  2 

(k-7o)J  -70)2> 


27o).o)2 


2  2 

270J.O32  (k- TWj-  TWg) 


2  toj-o^ 


2  To)^o>2 


■O). 


(k-Tdlj-Tft||) 


(A-82) 


Let  the  denominator  of  the  determinant  be  known  as    A  or  (A-82) . 


A  =  (k-70)j-  7o)2) 


2  2 

(k- 70)J  -  7o)2)        -270)10)2 

2  2 

27o)1o>2  (k- Tct)1  -  7o)2) 


-0)» 


2  2 

(k-  70)^  -  70)2) 


+   0), 


"2 

-270)^0)2 

"wl 

wl 

(k-  70)2-  70)2) 

-0)2 

2  TuijOij 

w2 

(k- 

2           2 

70)j  -  70)g) 

<"2 

(k-  70)2-  70)g) 

-Wj 

«, 

2  7o>.o>2 

"w2 

2  7O).0)2 

wl 

(k- 

7  0)2  -  7  0)2) 

2  70)-o)2 


0)r 


w. 


2  70)^0)2 


2  2 

(k-  70)J-70)2)       -27  0)1o)2 

2  2 

270)1O>2  (k-TWj-TWg) 


1 


'2 

(A-82) 
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2  2     I  2  22  2  2  22 

A  =  (k-  tw.  -  too%)     (k-Tco..  -  tco2)    +  w^k-rw.-TWn)  -2  rco.  co2 

+  o)2(k-  rco.  -  to>2)  +  2TW.  a)    +  4r    co..  co2(k-  rw,  -  rco2) 

[/■  2  2 \  ^  3        r>        2        /,  2  2x 

co2(k-  ro).  -  TWn)     +  ^o  +  2  r  co.  co2(k- r  co.  -  rco2) 

.223  2  0        2       „  2  2X1 

+  4r    co.  o)2  -  a).  a)2  +  2  r  co.  co2(k-  rco.  -  rco2)  I 

[2  2  2  2  2  2 

2  rco.co2(k-  rco-  -tco2)  +  co.  co2  -  co.  (k-  r  co.  -  r  co2) 

-  2  tco.  a>2  (k-  rco.  -  tco2)   -  co.  +  4r    co.  co2   I 

[3  2  2  2  2 

2toj1o>2  -  co..co2(k-  ret)--  tw2)  -  co.co2(k-  r  coj  -  t  co2) 

2  1 

+  2ra>1co2(k-  tco*-  t  a|)      -  2t  Wjcu2  +  8T3cOjCoi*  (A-82) 

A  0 

A    =  (k-T  Wj-TU>2)      +  (2a)j  +  2a)2)(k-Tu)j-T  w2) 

+  8r^co2(k-  TWj-T^)   -  2^w2  +  co*  +  a>2  -  16  f?C0j  co*         (A-82a) 

Determinant  Form  of  Coefficients  D,  E,  F,  and  G. 

0  -a>2  -co.  2rco.co2 

2  2 

0  (k- rco.  -  rco2)       -2x0^0*2  -co. 

2  2 

Aacoj  2  7u>.co2  (k-  r  coj  -  t  co2)  -co2 


D  = 


2  2 

Baco..  co.  co2  (k-T  ojj  -  r  co2) 


(A-81) 


E  = 


2            2 
(k-  T  CO-  -  r  co2) 

0 

co2 

0 

cox 

Aaco. 

2  TCO.COg 

Baco. 

-co.  2rco.co2 

-2tC0.C0«>  -co. 

2  2 

(k- t  co.  -  r  co2)  -co2 

2  2 

COg  (k- T  COj  -  r  co2) 


(A-83) 


82 


2            2 

(k-  liij-  Tiiij) 

-w2 

0 

"2 

(k-    TWj-T  Wg) 

0 

0)j 

2  TCJjO)2 

Aao) 

2  tcj^^2 

wl 

Bao> 

2  Taj.  c^2 


-Ci>< 


r.  2  2v 

(k-   TO)j  -  TO^J 


(A-84) 


2  2 

(k-  t  Wj  -  t  o>2) 


■o>* 


•  £J 


u>. 


2TWJOJ2 


(k-TWj-TWg)  -2TWJCJ2 

2  2 

2     T0)1Ci>2  (k-TCJj-  tu2) 


W1 


c>< 


0 

0 
Aaw, 
Baw. 


(A-85) 


AD=  Aao^   pTO^c^k-Tc^-Twl*)  +  4r    c^o^  Wj-WjWg 

+  2r  WjWgCk-  TWi  -Tw2)  +  a)1(k-Ta)1-rw2)  J 

2 
-Bau^   [-2  7  0)^)2+  210)^2+  21^(1^-1  u>x  -  tu>2) 


(A-81) 


900         3  2 

+  4t    OJjWg  +  <^i  ~^iw2 


[2  2  2 

4TW1aJ2(k-T  Wj  -t  a>2) 

+  Wl(k-TwJ-Tw|)  +    Baw1|2T«1w|-2TwJw2 

-  2  t  WjWg  (k-Ta)|-TW2)     -  8  T  wj  Wg  +  2wiC4,2(k  "  T  wl  "  T  W2^J 

(A-81) 
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[223  3 

-2TCJ1a)2(k-  r  ctJ-  -  r  w2)  +  2  r  w-cjg  +  2  tw.ho* 

q    o     o  2  2  2  2 

+  8r  Wj,  o>2  +  o>- Wofk-  tojj  -  tw«)  +  w^w^Ck-  to*..  -  tw^) 

[2  2  2  3  2 

2  toj- aj„  (k  -  tw.  -  to)2)  +  0)-  +  2  TWjOJg  (k 

232  2  2  2  1 

+  4r    WjCdg  +  cu^k- rwj  -  tcj2)      +  ^1^2  (A-83) 


2  2* 


3m3 

w2 


[22  3  3  *}     3 

2rco1a>2(k-  r^J  -  rwp     -2  tcj^2  -  2  TWjO;^  -  8  t°o>J 

22"1  T  2  223 

-2o)1o;2(k- to).  -  rcjg))      +     Baw*  I  4ro)1a;|(k-  rwf  -  T««)  +  wi 

232  2  2  21 

+  4t    oj.  0^2  +  cjjCk- rcjj  -  tcu2)      +  WiW,  (A-83) 


2    2 

>lw2 


[2  2  2     2 

(k-  tWj  -  t«2)      +  2  tcj*  a>2  +  2 

4  La  Li  Li    f  t  Li  Li   \  Lj    /  ,  Ld  Li   \  Li    f  i  A  Z   \ 

T     0)^  C4>2  (k-  T0>j  -  T0>2/  +  CJj  (k-  TWj  -  TOJg)  +  ^2^k~  ra,l  ~  TCt,2' 

r  2  22  223 

-  Baw..      -a>2(k- tciJj  -  rw2)     +  Wj  u>2  +  4  t    ^1^2 

2  2  2  2  2  2  3l 

-  2  TWj  WgCk-  TO),  -  TW2)  +  2  TCJ1  Ct>2(k-  TWj  -  tcj2)    -  k>2 


2.2 


2.2.2 


AF=     Aaa,1[(k-r^-TaJ2)     +  4  rc^  o>2  +  (o>x  -  4  r    c^  o>2  +  o>2) 


(A-84) 
.2- 


(k-  rojj  -  tco2) 

.223        3 

4  r    OJj  o>2  +  Wg 


[2  2    ■         2 

(^(k-  tw.  -  TOJg)      -  cj.  cog 


(A-84) 


AG  =  -Aaw 


4  r2  o>2  w3 


wlw2 


r  22 

.     a>2(k-  rcjj  -  to>2) 

2  2  2  2  2  2  3~l 

i-  2  rojjWgik-  rcjj  -  t<a>2)  +  2  rcuj  o>2(k-  tu^  -  tw2)  +  o>2 

t2  2  22  222  2  2 

(k  -  rojj  -  T<x>2 )      +  2  T(j)^  k>2  "  2  TCJi  ^2  +  wl  ^k  ~  ra,1  "  Tw2 ' 

>2-ra>2)]  (A-85) 


4  Li  Li  L*    ,  1  O  U  \ 

T     Wj  0>2(k-  TOJ^  -  TC4J2) 


o>2(k 
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[2  2  2  2     2 


4  TO) ^  OJgCk-  TO).  -  toj! 


I)  -«|1    +    Bawjhk-TwJ-TWg) 


+  (ojj  +  4  r    Wj  o>2  +  ^2  ^k  ~  T  wl  "  T  ^2  M 

(k  -  aw,  cos  w-t)  x.  +  x,  +  r  x.  =  (aw,  cos  o>,  t)  x« 


Let 


x5  +  x6 


x«  =   D  sin  w,t  sin  wyt  +  E  sin  w,t  cos  out 
+  F  cos  w,t  sin  wit  +  G  cos  w-t  cos  w2t 


where: 


(A-85) 

(A-66a) 
(A-86) 

(A-75) 
(A- 63) 


kx5  +  Xc  +  rx5  =  (aw,  cos  w,t)  x« 
Substitute  (A-86)  and  (A-75)  into  (A-66a). 

(k  -  aw,  cos  w,t)  (x&  +  xJ  +  (x&  +  xA  +  r (x&  +  xA 

=    Daw,  sin  w,t  cos  w,t  sin  wA  +  Eaw,  sin  w,t  cos  w,t  cos  wA 

2  2 

+  Faw,  cos    w,t  sin  wA  +  Gaw,  cos    w,t  cos  wA  (A-87) 


Try  a  solution  for  xX  where: 
Daw, 


Eaw. 


kx5  +  x5  +  rx5 


sin  2w-t  sin  Wgt  + 


sin  2a),  t  cos  w0t 


Faw. 
2 

Gaw, 


sin  (j)A  + 


cos  wA  + 


Faw, 
2 

Gaw 


cos  2w,t  sin  w2t 


-  cos  2w,t  cos  <j)A 


(A-88) 


Let 


=   H,  sin  2w,t  sin  wA  +  I  sin  2w-t  cos  uA 
+   J  cos  2w,t  sin  w2t    +  L  cos  2w,t  cos  uA 
+    M  sin  w2t  +  N  cos  wgt 


(A-89) 
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Xg  could  be  written  in  an  equivalent  form  with  certain  terms  containing 
(2o>1t  -  w2t)  by  use  of  trigonometry  relations.    Thus  a  time  dependent 
drift  term  is  seen  to  appear  when  the  harmonic  situation  2a)«  =  a>2 
exists. 

Xc  =    2H.,  <j)«  cos2w1tsinw2t   +    H..  u>2  sin2o>..tcos  aj„t 
+  21  Wj  cos2wjtcosa)2t     -    Ico2  sin2w1t  sinaj2t 

-  2Jcu^  sin2co1tsina)2t     +    Ja>2  cos2a>]tcos  a>2t 

-  2Lc^  sin 2a). t cos  uJl    -  L  a>2  cos 2 cousin a>2t 

+    Mw2cosa)2t    -    N  a>2  sin  ct>2t  (A-89a) 

Combining  terms  of  x&,  one  obtains: 

Xr  =  (2H^o)j  -  La,o^  cos  ^wi *  sin  w2*  +  ^1^2  "  2La,i)  sin  2w jt  cos  a>2t 
+  (2Icu-|  +  JwJ  cos2a>,tcosoJ2t    -  (Ia>2  +  2Ja>.)  sin2u>jtsin  w2t 
+  (Mcu2)  cos  w2t    -  (Na>2)  sin  a>2t  (A-89a) 

x5  =  (2H.0;.  -  LoO  (-2aO  sin 2a;. t  sin  o>2t  +  (2H-.0).  -  Lci>2)  (ct>2)  cos  2u)A cos  a>2t 

+  (H-.ct)2  -  2LaO  (2aO  cos2a)1tcosa>2t    +   (H.a^  -  2Lu>..)(-ci>2)  sin2aj«t  sina^t 

+  (2Ia>..  +  Ja>2)(-2w1)  sin2cJ1tcosct)2t   +  (21a).  +  Ja)2)(-a)2)  cos  2w1tsinw2t 

-  (Ia>2  +  2JuO(2aO  cx^oj-tsinajnt    -  (Ia>2  +  2Jct>1)(a>2)  sin2a>..tcos  a^t 

-  (Ma)2)sinw2t    -(Nw2)  coscJ2t  (A-89b) 

Combining  terms  of  x5,  one  obtains: 

x5  =  (2Lctf..w2  -  4H..  Wj  -  H.  u>2  +  2La>1a>2)  sin  2cu.,t  sin  a>2t 

2  2 

+  (2H1w1a)2 -Lct>2  +  21^(^0)2   -  4L(jj^)  cos  2ct>1t  cos  a>2t 

2  2 

+  (-4IW.    -  2Jw1o)2  -Iw2  -2Jo.)1ct)2)  sin  2a>jt  cos  o)2t 

2  2 

+  (-2Ia>.,u>2  -  Ja>2  -  2Ia)1a;2  -  4Ja>j)  cos  2w.t  sin  w^t 

2  2 

+  (-Mw2)  sina>2t   +    (-N«2)  cos  a>2t  (A-89b) 
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Now  substitute  (A-89),  (A-89a),  and  (A-89b)  into  (A-88)  and  solve  for 
coefficients  H^I,  J,  L,  M,  and  N. 

kH^  sin2w1tsinc4)2t  +  kl  siii2uA cos a>2t  +  kJcos2a><t  sina>2t 
+  kLcos2ct).jtcosa>2t  +  kMsincjgt  +  kNcosa^t]  + 
foHi^j -LaJ2)  eos2a)1tsina.)2t  +  (H..a>2  -  2LO)-)  sin2oJ1tcosa)2t 
+  (21  d)^  +  Ja>2)  cos  2a>jt  cos  a>2t  -  (I  a>2  +  2J  aO  sin 2cjA  sin  a>2t 
+  (M  a>2)  cos  Wgt  -  (N  cj2)  sin  w2t  J    + 

r  22 

r  (2L  0*^2  ~  ^1  wi  "^1^2  +  2La^Ct)2)  sin2aJ1tsina)2t 

2  2 

+  t(2H.  WjWg  -  La»2  +  2H.  aKa>2  -  4L  oj-  )  cos  2a>..t  cos  a>2t 

2  2 

+  t(  -4Iw-  -  2Jfc>1a>2  -  Ia>2  -  2Ja>,Ci>2)  sin  2a)1tcosa)2t 

2  2 

+  r(  -2Ia>.a>2  -  Ja?2  -  2IcjjO)2  -  4JuO  cos  2cjjtsin a>2t 

2  2 

+  t  (-M<x>2)  sin  cu2t  +    t(-Ncj2)  cosc^tl 

/  Daan \  /  Eaw.\ 

-I  sin  2aj|t  sin  a>2t  +  I  1   sin  2u> 


\    2 

Fao) 


jt  cos  out 


(Gaw. 


sin  ojgt   + 


COS  CJgt    + 


Faw 


Gaa). 


—J  cos  2a)jt  sin  a). 


cos  2co..t  cos  w2t 


(A -90) 


Now  equate  like  terms  of  the  various  trigonometric  functions  and  solve 
for  Hp  I5  J,  L,  M,  and  N.     (Note:   As  B,  C?  D5  E,  and  F  are  known) 


sin  2uA  sin  u)2t 


2  2 

kflL  -  (Io>2  +  2J  <*>,)  +   t(4Lo)1o;2 -4HJOJJ  -HjOJg    =    - 

2  2 

H1(k-4ru;1  -tcj2)  +  (-a>2)I  +  (-20^)  J  +  ^rc^a^)  L 


Daw, 


Daw, 


(A-91) 
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sin  2a)  A  cos  qoJl 

9  9  Eact)., 

kl  +  (H1o;2-2Lco1)  +    t(  -  4Ia>j  -  4^0^ -Ia>2)    - 


2 
H1  (oj2)  +  I(k-4rWj  -rcjp  +  J(~4TCt>1w2)  +  L(-2w-)    = 


2  2  EaWl 


2 

(A- 92) 


cos  2o)A  sin  a>2t 


Fawl 


kJ  +  (2Hu>..  -  Lai.)  +  t(-  4Io>1a)2  -  Ja>2  -  4Ja>.. )    = 

z 
9  9  Fao). 

H(2w1)  +  I(-4Tw1w2)  +  J(k-rct>2  -  4T  o>j)  +  L(-o>2J 


cos  2a>*t  cos  a>2t 


2 

(A- 93) 


2   .  0„..  ..        .T  ,.2,         Gaa,l 


2 


kL  +  (21a;..  +  JcuJ  +    r(2Ha)1C4)2  -Lw,+  2H  cUjO^g  -  4L  wf) 

9  9  Gaa?1 

H^toj-wJ  +  1(2^)  +  J(w0)  +  L(k-TO)o  -4rw7)    = 

11  l  l  *  1  2  (A-94) 


cos  u)Jl 


9         Gaw..  9  Gao)1 

Mw2-tNco2    =    L  M(w2)  +  N(-t«2) 


sin  a>2t 


2 
(A-95) 

Fact?!  0  Faw- 

*  _  -i  Tk/T  i     ~ ,.^\   .   XT/   ...  \    _  -i 


Noj2-tMo)2=    i-  M(-TWg)  +  N(-w2) 


2 
(A-96) 


Solve  (A-95)  and  (A-96)  for  M  and  N. 


Gao)., 

M    =    rNo)0  + i.  (A-95a) 

2         2^2 
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Substitute  (A-95a)  into  (A-96). 


-  Na>o   -  T  Qi)\ 


tNo^2  + 


Gaw. 
2o>0 


Fao). 


2 

Faw 


9    «       Graci).,^ 

Nc^-Nr2^3- -L±     '    -J 


N(-o9-t2w3)    = 


Faa>«  +  rGaoj-o) 


1~2 


N 


aw^F  +  tGcj2) 
-2w2(l  +  r2w2) 


N   = 


-aa^F  +  tGojJ 
2cj2(1  +  r    a>2) 


(A- 97) 


N    =    -t  Ma>« 


Faw. 
2a>0 


(A-96a) 


Substitute  (A-96a)  into  (A-95). 


Mw2  -  tcj2 


-tMoj0  - 


Fao). 


2wc 


Gacu. 


2M    3       TFawiw2 
Mu>2  +  r    Ma>2  +  


Gaw. 


(A-95) 


M(ct>2  +  r     a>2 


Gaw.  -  rFao^cjg 


M 


aco.  (G  -  r  Fo>2) 
2cj2(1  +  t    a>2 ) 


(A-98) 
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(k-4TWj  -  tu2)H1  +  (-uj2)I  +  (-20)^  +  (4rci)1oj2)L 
(w2)H1  +  (k-4ro^  -  tw^U  (-4tw1oj2)J+  (-20^)  L 


Dacj. 


Eao), 


9             9                               Fao)! 
(2w1)H1  +  (-4  70)^2)1+  (k-Tc^2  -4twJ)J+ (-oj2)L    =   i 


9            9              Gaoj, 
(4rc4)1w2)H1  +  (2^)1  +  (o?2)J+  (k-TW2-4TO)j)L    -  L 


(A-91) 


(A-92) 


(A-93) 


(A- 94) 


Solve  for  Hp  X,  J,  and  L  by  determinants. 
Daco., 


Hr 


2 
Eao? 


(-«,) 


(-2^) 


2 

Fao) 


1  2  2 

-    (k-4rct>1       -  tcj2)         (^tokwJ 


(A- 99) 


(-2^) 


Gaco. 


2  2 

(-47-0^0*2)  (k-4ro)1       -tw2)  (-w2^ 


(2^) 


(c2) 


2  2 

(k-  470^  -  7W2 ) 


2  2 

(cj2)  (k-47o?1  -  7w2 )       (-470)20)2) 


(-20)^ 


(470)^2) 
(-20^) 


(2o)x) 


2  2 

(-47a)1o)2)        (k-470).  -  7cj2)  (-^2) 


(2wx) 


2  2 

(u)2)  (k-470)1  -  70)^) 

(A- 100) 
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Let  the  denominator  of  the  determinant  be  called  X  or  (A- 100). 


1  = 


L  = 


(k-  4To)j-T 

«2> 

Dawj 
2 

(-20^) 

^TWjWg) 

(o>2) 

Eaw, 
2 

( 

-bru^o^) 

(-2&IJ) 

tH 
O 

1-H 

(2o)x) 

Fao>.. 
2 

(k- 

/1       2           2, 
-47o,1-70>2) 

(-c2) 

< 

{At(x)1U)2) 

Gao). 
2 

(cu2) 

(k 

2          2 

-47o>1  -  Two^ 

X 

2            2 

(k-470)j  -  TU)*) 

<-«2> 

Daw.. 
2 

(4x0)^2) 

(o>2) 

(k- 

AT     2 

70,2) 

Eao^ 
2 

(-20^) 

CM 

0 

(2wj) 

( 

-  4toy 

*2> 

Fao)j 
2 

(-u2) 

1 
< 

(47  w^) 

(2^) 

Gaojj 
2 

(k- 

^    2        2* 

-  47o>.  -  OJg) 

X 

0 

(k-4ro)1-r 

»|) 

<-«2> 

(-20^) 

Daojj 
2 

(a>2) 

(k- 

4TWJ- 

T«g) 

(-470)1oj2) 

Eaojj 
2 

O 

H 

(2Wl) 

( 

-470)jOJ 

2) 

(k-4ra,2_7o)2) 

Faojj 
2 

< 

(470^0*2) 

(20)^ 

<«2> 

Gao), 
2 
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\=  (k-47Ci)2-70>2) 


(k-  4to)2  -7o>2)        (-4tWi«2) 


(-470^0)2) 


(-«,) 


(2^) 


(o>2) 
(20^) 

(470)^2) 


(w,) 


+  (-2^) 

(2^) 

(4  7Wlw2) 

(0>2) 

470)^2) 

pwj) 

^TWjWg) 

(k-  47o)j  -  T^2) 
(«2> 


(-470)^2) 

2  2 

(k-47cu1  -  7o>2) 

(Wo) 


2  2 

(k-47cu1  -  7Ct)2) 

(-470)^2) 

(2CDJ 


(-2^) 

(-o,2) 
(k-4TWJ-Ta)2) 


(-2^) 

(-o>2) 

2  2 

(k-47w1  -  tw2) 


(-o>2) 

2  2 

(k-47o>1  _  7^2) 


2  2 

(k-4r<xJ1  -  70>2)        (-470)^2) 


(2^) 


(k-4  70)1  -  TO),) 


1 

w 


(A- 100) 


X  =  (k-47a>2-7a)2)  [(k-47a>2-7a|)     +  (Stu>\ w|)  +  (8ra£a|) 

-(-4a>j)(k-4TO)j-Tw|)  -  (16  72  cuj  a)2)  (k- 470)2 -70)3) 

-(-w2)(k-47cJj[-7a,2)] 

2 
+  (w2)  [(co2)(k-47U)2_Ta)2)     +  (1672a)2^)  -  (4a)2a)2) 

-  (-87a)2a)2)(k-47a)2_Ta)2)  -  (-87a>2CD2)(k-47a>2  -  tw\)  +  (a>|)] 

r  2  22  222 

+  (-2a),)     (-47W1a)2)  (k-470)1  -r^)  +  (-470)1a)2)(k-  4t/«7  -  7Ct>2) 
+  (-8wJ)  -  (3272a)Ja;2)  -  (2Wl)(k-  4Ta)2  -  7a>2)2  -  (-2a>1a)|)] 

r  3  22  3 

+  (-47a)1w2)      (-470)^2)  +  (4Ta>1a)2)(k-4Ta)1 -7  w2)     +  (-16twJ  u>2) 

-  (64  73o)Ja)^)  -(2a)1w2)(k-47a;2-7a)2)-(2a)1W2)(k-47a)2-7a)2)] 

(A- 100) 
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r  4 

=    F(k-T«*-Ta|)     +  (8o>2  -32  r2  o>2  o|  +  2o>2)(k- 47o>2  -  7o>2) 

+  (64  7o>2o>2)(k-47o>2-7o>2)  +  (o>*  +  16o>*  +  32r2w2^ 


+  128  72o>*o>2  +  256 u>*w*  T4 


-8c2^2J 


(A- 100) 


Daw. 


AH. 


2  2 

(k-4rcj1  -  ro>2)        (-470)^2) 

2  2 

(-470^0*2)  (k-47o>1  -  7o)2) 


(2^) 


(«2) 


(-20)^ 

(-<*2) 

2  2 

(k-47w^-7a)2) 


Eacj, 


(-c2) 

(-470)^2) 
(2^) 


(-2Wl) 


k-47o> 


7a|) 


(470^^2) 


(o,2) 


(-co2) 
2 


(k-47co 


ro-2) 


Fao>. 


(-cu2) 


(2Wl) 


(-20^) 


2  2 

(k-47^1  -  7w2)        (-47  0)^2) 


K) 


(47  o)^) 

(-20^) 

2  2 

(k-  47  w.  -  7oO 


Gao> 


(-o>2)  (-20^) 

2  2 

(k -  47  0).  -  70)2)        (-470). o^) 

2  2 

(-470)^2)  (k-4ro)1 -70>2) 


(470)1o)2) 

(-20^) 

(-o)2) 

(A-99) 
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Hx  =  I -lj    J    (k-4TwJ-Td|)      +  {%T(ti\u\)  +  (8to>2o|) 

+  (4a>2)(k  -  4to>2  -  tu>\)  -  (16r2  u>2  w|)(k  -  4ro>2  -  Tu>2J 

+  (ct>2)(k-4rw2-Ta;2)l 

|_awij       (w2)(k_4rw2^rw2)     _(4a,2W2)+  (16t2w2w3j 

2  22  2  2231 

+  (8TW10J2)(k-  4x0^  -  toj2)  +  (Stco.Wq)^-4'''^!  -  TCt>2^  +  ^w2^ 


>2)(k-4Ta>2-Ta>2)  +  (8aiJ) 


+  (4rw1a)2)(k-4r^-TCJ2)  +  (32r2  ojJ  oj2)  +  (2a)]L)(k-4rw2  -  rw2) 


-  (2wlW|)] 

+  (64t3oj^^)  +  (4<w1w2)(k-4TWj-  r  a>2)J 


3  3  2  2 

«)  +  (1Qtu)^u)2)  -(4rcj1a>2)(k-4ra;1  -  tcj2) 


(A- 99) 
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*H,  = 


Daw. 


(k-4rw?-  two)     +  (16tcj?w2)  + 


1~2 


Eaw, 


Faw. 


+  (4ci>i-16t2wJci^  +  OJ2)(k-4r^-TW2) 


(a>2)(k-4T^-TCJ2)     -(4oj^2)  +  (16t2u>2u>;*)  + 


+  (16To,2o;2)(k-4T^-Taj2)  +  (a>i*) 


(2a)1)(k-4TCJ?-Ta;2)    -  (20),^)  +  32t2o>?o>2  + 


lw2 


1~2 


Gaa^ 


+  (8TaJ1a)2)(k-4ra>2-rct)2)  +  (8a>J) 


2  2  3 

(-4Ta>jO>2)(k-  4to>,  -  toO      +  4ra)..a>2  + 

3     3     3  2  2 

+  64r    a>j  w2  +  (4Tw1w2)(k-4Tw1  -  TWg)  + 


+  16tci)j  c^2 


(A-99) 
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Daw. 


AI=-| 


EaaK 


Fao). 


Gaco, 


(o>2)             ( 

-470)^2) 

(-2w2) 

(2^)        (k- 

2            2 

470J^-  70J2> 

(-o)2) 

(4tcj^oj2) 

(o>2) 

(k- 

2           2 
-470JJ  -  7w2) 

(k-4T^-T^) 

(-2^) 

(47u)1o)2) 

(2Wl) 

9 

(k-  47  Cl^  -  TU,) 

(-«2) 

(4ro)1o)2) 

(u>2> 

2          2 
(k-  47  ojj  -  T^g) 

2           2 

(k-47  0).  -  TW„) 

(-2^) 

(470)^2) 

(a>2> 

(-470)^2) 

(-2Wl) 

(47  w^) 

(c2) 

2           2 

(k-  TWj  -  7Ci>2> 

(k-  47  0)^  -  r^) 

(-2*^) 

(47(^^2) 

(o>2) 

(-470^0*2) 

(-2^) 

(20)^ 

2 
(k-  47  oj.  -  7C1O 

(-o>2> 

(A-101) 
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Dao), 


XI 


Eao>. 


Fao). 


(w2)(k-4TWj-TW2>    +  (16t2o)j^)  -(4o>?Ci>2)  + 

2  2  2  2  2  2 
+  (870^  w2)  (k-  4TOJJ  -  toj2)  +  (8raJ1a)2)(k-4Ta)1-  to>2) 

31 
+  (a>2)    + 

3 
(k-4ra^-Ta>2)    +(8r^a)2)  +  (8r^a>2)  + 

-  (16r2^^)(k-4r^-TaJ2)  +  (4o)2)(k-47  0)2-7  0)2)  -t 

+  (o|)(k-47  0)2-7O)2) 

2  2  3 

(-4Ta)1cu2)(k-4TCtJ1  -  tci>2)    +(16toj^  w2)  + 

3  3     3     3  2  2 

+  (47o)..o)2)  -  (64t    oj.  u>2)  +  (2w1a;2)(k-4TaJ1  -  70O  + 


rGao). 


2  2 

+  (2a)1o)2)(k-4rajj  -  tco2) 


(4t  a^  )(k-4Ta>j-Ta>2>  +  (8a)j)  + 

+  Hra^Mk^To^-Tal)  +  (3272o)Jo)2)-  (20)^^)  + 


+  (2o)1)(k-4r  Wj-rw^) 
(A-101) 
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\J   = 


/d 


ao). 


(a>2) 


2  2 


■Eao). 


2  2 

(k-4ro)1-  tcj2) 

(2wx) 

(47  0)^2) 


(-o,2) 

(-4x0^0*2) 

(2^) 


(4Tb>£U>2) 

(-«2) 

2  2 

(k-4roj1  -  rojo) 


Fao>. 


(k-4rCt)'1 
(«2) 
(4to)1o)2) 


™2> 


(-a,2) 

2  2 

(k-  4T0JJ  -  TOJg) 

(2^) 


(470)^2) 
(-2«x) 

(k-4TwJ-TW2> 


•Gaw, 


(a>2> 
(2WJ) 


(-a>2> 
2 
1" 


(k-4ra>i-To^)      (-2^) 


(4rw1w2) 
[-2&IJ) 


(A-102) 


Daw. 


XL 


(cu2) 


2  2 

(k -  rcjj  -to),)        (-4x0)^2) 

2  2 

(-470^0*2)  (k-4rw1  -  TCJ2^ 


(2wJ 


(w2) 


Eao). 


(k-4TWj-TW*) 

(2^) 


(-oi2) 

(-470*^2) 
(2W1) 


(-2^) 

2  2 

(k-470).  -  7Ct)2) 


Faco. 


(k-470^  -  7w2) 

(co2) 


(-«2> 

2  2 

(k-47^1 -70>2)     (-47a)1cu2) 


(2^) 


(-2^) 

t70)j( 

(«2) 


Gaa>. 


2  2 

(k-  47CJ..  -7Ct)2) 

(u>2) 
(2wx) 


(-w2)  (-2wj) 

2  2 

(k-  470),  -  7CJ2)        (-470^0*2) 

2  2 

(-47CJ.0J2)  (k-4Tw-  -  7a>2) 


(A-103) 
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A  J  and  AL  are  left  in  the  form  of  a  third  order  determinant.    Fur- 
ther determinant  breakdown  will  yield  terms  similar  to  those  associated 
with  AH  and  A  I. 

The  determinant  form  does  permit  the  detection  of  the  common 
factor  "a"  in  all  the  coefficients.    As  in  the  case  for  zero  elastic  re- 
straint, terms  beyond  Xr  may  be  neglected  if  the  order  of  magnitude 
of  "a"  is  much  smaller  than  unity.    Each  successive  term  in  the  series 
of  x  is  smaller  than  its  predecessor  in  accordance  with  the  value  of 
"a". 

x  =  x1  +  x3  +  x5+.  .  .   .  +x2m_1  +  x2m  (A-3) 

x-  =  A  sin  wJl  +  B  cos  a>2t  (A- 67) 

x«  =  D  sin  u)A  sin  o^t  +  E  sin  u)A  cos  ct)„t 

+  F  cos  a>-.t  sin  u<A  +  G  cos  wA  cos  wJl  (A-75) 

Xr  =  HjSin  2(j)A.  sin  a>2t  +  I  sin  2ooA.  cos  a>2t 
+  J  cos  2a>jt  sin  ct>2t  +  L  cos  2a).t  cos  a>2t 
+  M  sin  a>2t   +    N  cos  a>2t  (A-89) 

Drift  terms  may  be  shown  to  exist  by  applying  trigonometry  re- 
lations to  x«  and  x,..    Equal  input  frequencies  and  harmonics  will 
cause  time  independent  terms  of  constant  angular  displacement  to 
appear  but  there  are  no  time  dependent  terms  as  in  the  case  of  zero 
elastic  restraint. 
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x3    =    - 
6         2 

2 

F 

+  — 

2 
G 


C0S(0)jt-  Ci^t)   -  COS(cJ..t  +   CJ«t)J 

sin(c«;^t+  Wgt)  +  sinCcu.t  -w^t)] 
sin(o>jt  +  o)2t)  -  sinfw-t-Wgt)  ] 
cos  (a>- 1  +  Wot)  +  cos  (w.  t  -  wJl)  J 

cos  (2w1t  -  o)2t)  -cos(2ct>1t+  c^2t)  J 

sin(2u>jt  +  co2t)  +  sin(2w1t  -k>2t)  ] 

sin(2a)-t  +  o^t)  -  sin(2ct);.t- o>2t)  ] 

cos  (2a>]t  +  (j)  J)  +  cos  (20)-!  -  ajgt)  ] 


(A -7  5a) 


+    M  sin  u)Jl  +    N  cos  u)Jl 


(A-89a) 


x.  terms  are  pure  sinusoids  and  the  time  average  value  at  the  end 
of  complete  cycles  is  zero.    x«  terms  are  sinusoids  which  do  not 
display  time  dependent  characteristics  if  cu*  =  u>2.    Constant  dis- 
placements such  as  D/2  and  G/2  do  appear  when  u)«  =  a>2.    Likewise 
in  Xc  terms,  time  dependent  terms  are  lacking  when  2u>..  =  u>2  but 
constant  angular  displacements  such  as  H/2  and  L/2  are  present. 
This  phenomena  can  be  easily  explained  through  the  coefficients 
attached  to  the  sinusoids.     The  denominators  of  such  coefficients 
become  zero  for  w*  =  oj„  an<^  2w1  =  a>2  in  the  case  for  zero  elastic 
restraint.    However,  a  finite  elastic  restraint  keeps  a  real  value 
attached  to  the  coefficients  and  thus  indeterminate  forms  do  not  ap- 
pear when  sinCc^t  -  a>2t)  and  sin  (2^  -o>2t)  are  zero. 


101 


This  analysis  is  verified  by  Ref.   10  for  the  case  where  u*  =  cj« 
and  k  is  finite. 

One  other  possibility  is  seen  to  exist.    It  concerns  A  and  A , 
the  denominators  associated  with  the  determinant  forms  of  the  co- 
efficients.   If  A  and  A  can  be  made  equal  to  zero  by  some  combination 
of  q)a  and  Wg,  drifts  could  possibly  occur.    Attempts  to  solve  A  for 
(j)+  to  give  A  =  0,   lead  to  "an  equation  in  u)«  of  the  eighth  order.    Simi- 
lar difficulty  appears  when  A  is  set  equal  to  zero  and  an  attempt  is 
made  to  solve  for  u)..    With  such  values  of  a*,  and  a>„  which  satisfy 
A  =  0  or  A  =0,  L'  Hospital's  Rule  could  be  applied  and  the  possibility 
of  a  drift  rate  could  be  investigated.    This  appears  as  a  field  for  more 
complete  investigation. 
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